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Summary 

A full aero-thermal calibration of the NASA Glenn Icing Research Tunnel was completed in 2004 
following the replacement of the inlet guide vanes upstream of the tunnel drive system and improvement 
to the facility total temperature instrumentation. This calibration test provided data used to fully document 
the aero-thermal flow quality in the IRT test section and to construct calibration curves for the operation 
of the IRT. The 2004 test was also the first to use the 2-D RTD array, an improved total temperature 
calibration measurement platform. 


Introduction 

Updating and maintaining the icing cloud and aero-thermal calibration and flow quality data are 
critical tasks in the NASA Glenn Icing Research Tunnel (IRT). Establishing and maintaining the 
calibrations is the only way to insure that correct test section operating conditions are being set for a given 
set of facility control inputs. Without knowledge of the calibration and flow quality data, it is impossible 
to provide high-quality test data for customers and researchers. To insure that the data quality from the 
IRT is always at the highest standards, IRT staff and wind tunnel calibration personnel have developed 
and continue to refine the calibration methodologies used. 

The aero-thermal calibration test program provides the calibration relationships needed to accurately 
determine the test section operating conditions for total pressure, airspeed (comparable Mach number) 
and total temperature. Once these calibrated conditions are determined, all other aero-thermal operating 
conditions can be calculated. In addition to the calibration information, detailed flow quality data (spatial 
uniformity) is also collected. The collected calibration and flow quality data provide a complete 
description of the aero-thermal operation of the IRT. 

The current aero-thermal calibration methodology and instrumentation has evolved based on lessons 
learned by the calibration team over the past 15 years, input from IRT customers, in-house research staff, 
IRT facility personnel, and recommendations from References 1 and 2. Reference 1 defines a time table 
and level of calibration that should be conducted following modifications to an icing tunnel. 
Correspondingly, Reference 1 recommends that short check calibrations should be conducted multiple 
times each year to insure repeatable operation of the wind tunnel. Based on these recommendations, the 
following aero-thermal calibration tests have been defined: 

• Full calibration: Includes collection of aero-thermal calibration data over the entire operating 
range of the IRT as well as a complete mapping of the IRT flow quality at one cross section of the 
test section, typically at the centerline of the turntable. A full calibration should be conducted 
every 5 years or following a major facility modification. The full calibration provides all of the 
information needed to develop the aero-thermal calibration relationships, or calibration curves, 
for the facility and provides extensive detailed information on the aero-thermal flow quality 
characteristics of the test section. 
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• Interim calibration: Includes collection of aero-thermal calibration data over the entire operating 
range of the IRT, however, data is only required from the test section centerline. An interim 
calibration is conducted one year following a full calibration. The interim calibration is used to 
validate the data collected during the full calibration, therefore, if the full calibration and interim 
calibration data match then the facility is operating properly, however if not, then an investigation 
must commence immediately to rectify the discrepancies. 

• Check calibration: Ideally a one-day test using a very simple rake. Check calibration testing 
covers the entire operating range, but the test matrix is not as extensive as a full or interim 
calibration. This calibration should be conducted 2 to 3 times per year. The check calibration data 
is used in a statistical process control analysis to monitor the stability and repeatability of the 
facility operation. 

This paper provides a detailed description of the full aero-thermal calibration conducted in the IRT in 
2004. This calibration was required due to replacement of the drive system inlet guide vanes, which is 
detailed in the subsequent Description of Facility section. The 2004 full calibration was also needed to 
update the total temperature calibration based on improvements to facility instrumentation (Ref. 3). 
Previous to the 2004 calibration, the most recent full calibration was completed in 2000 following the 
replacement of the facility heat exchanger (Ref. 4). 

The goals for the 2004 aero-thermal full calibration are listed below: 

• Qualify for operation the 2-D RTD array as the new measurement platform for total temperature 
calibration and flow quality data in the IRT test section. 

• Complete data collection and construct aero-thermal calibration curves for the IRT. The three 
calibration relationships used are total pressure, Mach number and total temperature. 

• Complete detailed flow quality surveys using the 9-ft survey rake and 2-D RTD array. 

While the bulk of this report focuses on the full aero-thermal calibration conducted in 2004, a 
discussion of the 2005 interim calibration is included in an appendix. The 2005 interim calibration 
provided additional information on the total temperature calibration of the IRT and allowed an 
opportunity to verify additional planned improvements to the total temperature calibration. 1 

Symbols 


a 


Aq,A\,A 2 
Co, C u C 2 
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Ctt,o? Ctt,\ 
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K {) to K 2 
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P 


speed of sound, ft/sec 
Hot wire curve fit coefficients 

RTD probe total temperature flow recovery curve fit coefficients 
Pitch angle pressure coefficient 
Yaw angle pressure coefficient 

Test section total temperature calibration curve fit coefficients 
Total pressure coefficient (9-ft rake pressure probes) 

Static pressure coefficient (9-ft rake pressure probes) 

Total temperature recovery coefficient (9-ft rake total temperature probes) 

Hot wire anemometer output voltage, volts 

Flow angle prediction coefficients, degrees 

characteristic length, ft 

Mach number 

Number of data points 

Pressure, psia 


1 Check calibration data are collected as part of both full and interim calibration test entries, as well as during 
dedicated check calibration test entries. However, the check calibration tests are not covered in this report. 
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Flow angle probe pressures, psia 
Spray bar air pressure, psig 
Average of P u P 2 , P3, and P 4 , psia 
Ratio of static to total pressure 
Static pressure, psia (psf) 

Average of P 6 , P 7 , P%, and P 9 (probe static pressure), psia 
Total pressure, psia (psf) 

Equal to P 5 (probe total pressure), psia 
Dynamic pressure, psi 

Specific gas constant for air, 1716 ft 2 /(sec 2 °R) 

Reynolds number 

Root-sum-squared uncertainty estimate 
Temperature, °C (°F) 

Axial turbulence intensity, percent 
Transverse turbulence intensity, percent 
Static temperature, °C (°F) 

Total temperature, °C (°F) 

Average of the 24 D-comer total temperature measurements, °C (°F) 
D-comer total temperature measurement, °C (°F) 

Velocity, knots (ft/sec, mph) 

Hot wire axial velocity, ft/sec 
Hot wire transverse velocity, ft/sec 
Flow angle resultant vector magnitude, degrees 
Flow angle resultant vector orientation, degrees 

Axial coordinate with axis origin at bellmouth/test section weld seam, in. 
Spanwise coordinate with axis origin at the test section inner wall, in. 
Vertical coordinate with axis origin at the test section floor, in. 

Pitch flow angle, degrees 

Yaw flow angle, degrees 

Ratio of specific heats, 1 .4 

Measured pitch or yaw angle offset, degrees 

D-Comer average temperature, °C (°F) 

Standard deviation 

Standard deviation of the D-comer total temperature, °C (°F) 

Hot wire sensor inclination angle, degrees 

Parameter used to calculate flow angle resultant vector orientation, degrees 


Subscripts 


arm Pertaining to the 6 degree-of- freedom digitizing arm 

array Pertaining to the 2-D RTD array 

avg Average 

bm Bellmouth conditions 

bm, north Measured by the north bellmouth rake 

bm, south Measured by the south bellmouth rake 

eff Pertaining to hot wire effective velocities 

HW Pertaining to a hot wire probe 

i Data point index or hot wire sensor index 

level Pertaining to a digital level 
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local 


Pertaining to parameters in the test section after all probe calibration coefficients have been 
taken into account 
nom Nominal test section condition or setting 

rake Pertaining to the 9-ft survey rake 

ruler Pertaining to a tape measure or ruler 

ts Test section conditions (based on calibration data) 


Description of Facility 

The NASA Glenn Icing Research Tunnel (IRT) is a closed-loop atmospheric tunnel illustrated in 
Figure 1. The airflow is driven by a 25-ft diameter twelve blade fan that is powered by a 5000-hp electric 
motor. The test section, shown at the bottom center of Figure 1, is 6-ft high, 9-ft wide, and 20-ft long, 
with no divergence along any of the test section surfaces. The velocity in the empty test section can be 
varied between 25 and 340 knots at -17.8 °C (37 to 575 ft/sec at 0 °F). A pair of heat exchangers located 
in the C-D leg of the tunnel loop, seen at the left of Figure 1, allow for total temperature settings as low as 
-30 °C. A set of ten horizontally oriented spray bars, located in the settling chamber at the bellmouth inlet 
prior to the test section, inject atomized water into the airflow to create icing conditions. The IRT has 
been in service since 1944 and has undergone several major upgrades and modifications over the years to 
insure that it remains the lead ground test facility for icing research (Refs. 5 and 6). In 2004, the inlet 
guide vanes (IGVs) upstream of the drive fan were replaced, primarily for structural reasons, but also to 
improve aerodynamic performance. Figure 2 shows the final configuration of the replacement IGVs 
upstream of the IRT drive fan prior to being painted. The view in Figure 2 looks downstream from Corner 
B showing the inlet and outlet guide vanes, fan blades, and floor drain. During the downtime for the IGV 
replacement, the new facility instrumentation for total temperature measurement was installed in D- 
corner, consisting of 24 RTD probes. More information about the IRT is located in Reference 7. 



Figure 1. — Layout of the NASA Glenn Icing Research Tunnel. 
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Figure 2. — Installation photo of the replacement inlet guide vanes (IGV) 
upstream of the IRT drive fan. 


Test Hardware, Instrumentation and Data Systems 

The 9-ft survey rake and the 2-D RTD array are two pieces of test hardware used to complete the 
aero-thermal calibration of the IRT. The 9-ft survey rake is used to collect pressure, flow angle, 
turbulence calibration, and flow quality data. The 2-D RTD array is used to collect total temperature 
calibration and flow quality data 2 . Detailed descriptions of the hardware, instrumentation and data 
systems are included in this section. 


9-ft Survey Rake 

Figure 3 shows the 9-ft horizontal survey rake installed in the IRT test section at the vertical 
centerline, a height of Z=36-in. The flow angle pressure probes, total temperature probes, and hot wire 
probes mounted on the 9-ft survey rake in Figure 3, will be discussed in more detail later in this section. 
The reference axis is also defined in Figure 3 depicting the vertical direction as the Z-axis, the horizontal 
direction as the 7-direction, and the axial direction as the X-axis. The origin of the left-handed X-Y-Z 
coordinate system is as follows: the X axis origin is at the bellmouth/test section weld seam as seen in 
Figure 4, the 7 axis origin is at the inner wall as seen in Figure 4, and the Z axis origin is at the floor as 
seen in Figure 3. This rake is used to collect calibration and flow quality data for pressure, airspeed, flow 
angle and turbulence in the IRT test section. 

Figure 4 illustrates the plan view of the IRT test section with the survey rake installed. The rake was 
installed at an axial position of X=l 79.3-in. This was also the axial location of the cross-sectional plane 
surveyed during the 2004 test program as well as the test program in 1997. Figure 5 shows the plan and 
front views of the survey rake. The following illustration, Figure 6 details a cross-section of the survey 
rake with exploded views of the flow angle pressure probe ports. 


2 Prior to the 2004 calibration, the 9-ft survey rake was used to collect all calibration data, including total temperature 
data. Starting in 2004, an improved total temperature measurement platform, the 2-D RTD Array, was introduced 
for calibration and flow quality data collection. 
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Figure 4. — Plan view of the IRT test section showing the installed location of the 9-ft survey rake. 
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A-A B-B 

exploded exploded 

view view 

Total pressure port detail Static pressure port detail 

Figure 6. — Cross section of the 9-ft survey rake showing flow probes details. 

The instrumentation mounted on the 9-ft survey rake includes eleven flow angle pressure probes, 
eleven total temperature thermocouple probes, and three thermal anemometry probes. The spanwise, Y 
positions of all probes with respect to the test section inner wall are defined in Figure 5. The survey rake 
was supported at both ends by bolting the rake to wall mounted support plates as seen in Figure 3. These 
support plates are 6-ft tall and have a bolt pattern that allows the rake to be positioned vertically every 


NAS A/CR— 20 1 0-2 1 6733 


7 



6-in. above or below the test section centerline. A vertical stmt, shown in Figure 3, was also used to 
support the rake in the center. As seen in Figure 6, the main body of the rake is formed by two aluminum 
I-beams and thin aluminum plates riveted together. Thin sheets of aluminum were bent to form the 
leading and tailing edge of the rake body. 

The flow angle pressure probes are bolted to the rake body I-beams as illustrated in Figure 6. These 
probes are of hemi-spherical head design, bottom left Figure 6, and have five total pressure ports in the 
head. The bottom center Figure 6 shows the center port in the hemi-spherical head design which measures 
total pressure, and the four circumferential total pressure ports which measure pitch and yaw angle. Four 
static pressure ports are located downstream of the head depicted bottom right Figure 6. Top Figure 6, 
shows the heads of the probes located 20.75-in. upstream of the rake leading edge and the static pressure 
taps on each probe located 5. 75 -in. downstream of the head. The probes were calibrated for Mach 
numbers between 0.1 and 0.6 in the NASA Glenn 3.5-in. diameter free jet calibration facility. Details of 
the probe calibration can be found in Reference 8. 

Aspirated total temperature type T thermocouple probes with copper/constantan wires were used to 
monitor test section temperature. These probes are also depicted in Figure 3, Figure 5, and Figure 6. The 
probes are mounted to the bottom surface of the survey rake with the tips of the probes approximately 
0.5-in. upstream of the rake leading edge and 2.0-in. below the rake centerline. The total temperature 
probes were calibrated for total temperature recovery in the NASA Glenn 3.5-in. diameter free jet 
calibration facility. 

Three single sensor and dual sensor (cross-wire) hot-wire probes with 0.00015-in. diameter tungsten 
wires, were used to measure the turbulence levels in the test section. The probes were mounted to the 
upper surface of the rake body as shown in Figure 3, Figure 5, and Figure 6. The top of Figure 6 more 
clearly illustrates the probe tips were located 9.75 in. upstream of the rake body leading edge and 3.0-in. 
above the rake centerline. 


2-D RTD Array 

The 2-D RTD Array was designed and built to provide a more accurate total temperature calibration 
and temperature mapping of the test section, with the benefit of reducing overall test time and cost. The 
total temperature probes used on the array are 4-wire resistive temperature device (RTD) with a ceramic 
capsule sensor. Figure 7 shows the details of the RTD probe. See Reference 3 for a discussion of the 
thought process used to arrive at the design of the array and selection of the probes. 

Reference 1 was utilized to determine the number and arrangement of probes on the array. Figure 8 
shows the resulting layout is a 7-by-7 grid with the probes spaced at 0, ±25, ±50 and ±75 percent of the 
test section dimension from the horizontal and vertical centerlines. The intersections of the dotted lines in 
this figure indicate the probe positions. The 7-by-7 grid arrangement provided the test section coverage 
required without incurring an unreasonable cost for the RTD probes. 
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Figure 7. — Details of RTD probes used on the 2-D RTD array (dimensions in inches). 
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Figure 9 shows the 2-D RTD array installed in the IRT test section (views looking downstream and 
upstream). The structure of the 2-D RTD array includes a set of seven vertical members and horizontal 
members. Each vertical member houses seven RTD probes and instrumentation wiring, while the four 
horizontal members provide structural rigidity. The vertical members seen in Figure 9 are 0.75 in. thick 
by 4 in. wide and have a round leading edge. Figure 10 illustrates the cross section of an array vertical 
member assembly. Also, Figure 10 shows how the tips of the RTDs extend 7 in. forward of the leading 
edge, and the 0.75 in. deep channel that was attached to the trailing edge to house the instrumentation 
lines. The horizontal members seen in Figure 9 are 0.5 in. thick by 3 in. wide and also have a round 
leading edge. The array was secured to the test section walls, floor and ceiling, with allowances made to 
account for thermal expansion of the array. The tips of the RTD probes were positioned 194.5 in. 
downstream of the tunnel datum line or 121.5 in. from the inlet of the test section (the thermocouples on 
the 9-ft survey rake were 199.5 in. downstream of the datum line). 

36 . . . . . . . 


27 



-27 : -; : ; : ; 


-36 ! i i s : i 

-54 40.5 -27 -13.5 0 13.5 27 40.5 54 

Distance from centerline, in. 

Figure 8. — 7x7 grid arrangement of the RTD probes in the IRT 
test section. 



Figure 9. — RTD Array installed in the IRT test section, (a) View looking downstream, (b) View looking upstream. 
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Figure 10. — RTD array vertical member cross sectional details. 



Figure 1 1 . — Typical total temperature flow recovery curve 
for the 2-D RTD array probes. 



Mach number, M 

Figure 12. — Total temperature flow recovery data for 
RTD probes 1 through 23 with average curve fit. 


A portion of the RTD probes (serial numbers 1 through 23) were calibrated for temperature flow 
recovery in the 3.5-in. Diameter Free Jet Facility. The recovery testing was conducted to determine how 
well the probes sense the total temperature over a range of airflow conditions. During the testing, the free 
jet air temperature was set to a nominal condition of about 520 °R (60.33 °F) and held constant while the 
free jet Mach number was varied from 0 to 0.5. Figure 1 1 shows the total temperature recovery results for 
a typical probe (serial number 15) with its individual flow recovery curve. The recovery curves were 
subsequently used to correct data collected in the IRT. The 20 RTD probes that were individually 
calibrated for flow recovery were mounted in the center of the 2-D RTD array. The probes not measured 
directly for their recovery characteristics were mounted around the periphery of the array and used the 
average recovery from all the calibrated probes, seen in Figure 12. The probes have an absolute 
measurement uncertainty of ±0.10 °C (±0.18 °F) for a probe that was individually calibrated for total 
temperature flow recovery or ±0.16 °C (±0.29 °F) for a probe that was not specifically calibrated for 
recovery (uses the average recovery values). The details of this uncertainty analysis are found in 
Appendix A. A summary of the flow recovery correction coefficients is provided in Appendix B. 


Facility Instrumentation 

The following standard facility instrumentation was used during this test program. 

1) Bellmouth pitot-static probes: A probe is mounted on the south and north wall near the exit of the 
bellmouth. Transducers associated with these probes measure the total and differential pressures. 
These measurements are used to arrive at Ptm and P s ,bm • The location of these probes is shown in 
Figure 4, as the outer (north) probe and the inner (south) probe. 
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2) Total temperature probes in D-corner: Twenty-four RTD probes are arrayed on the leading edges 
of the turning vanes in D-corner to measure the temperature profile exiting the facility cooler. The 
distribution of the measured temperatures and the average of these twenty-four probes, T^ davg , is 
used to monitor the operation of the heat exchanger. The locations of these probes are shown in 
Figure 13. 


Steady-State Data Acquisition System 

Real-time steady-state data acquisition and data display are provided by the NASA Glenn Escort 
Alpha system. This is the standard data acquisition and data display system used in the large test facilities 
at NASA Glenn. The system accommodates inputs from the Electronically Scanned Pressure System 
(ESP), the facility distributed process control system, and any analog devices such as thermocouples, 
RTDs and pressure transducers. This system records all steady-state pressures and temperatures from the 
9-ft survey rake, 2-D RTD array, the tunnel bellmouth pitot-static probes, and the total temperature 
probes in the D-comer. It also records important facility parameters such as fan rotational speed and spray 
bar air pressure. The Escort program, designated D046, used for data collection, online data analysis and 
display was developed specifically for the aero-thermal calibration testing. 

The Electronically Scanned Pressure system used during this test program utilized plug-in modules 
that each contained 32 individual transducers with individual ports. Each transducer/port can be addressed 
separately and scanned at a rate of 10,000 ports per second. Calibration of all ESP transducers was 
performed automatically by the system at least every 2 hr. For this test program, ±5 psid modules were 
used. The accuracy in the pressure measurements made with the ESP system is 0.1 percent of full-scale of 
the module, or 0.005 psia for the ±5 psid modules used. 

Thermal Anemometry Instrumentation 

Commercially available constant temperature hot wire anemometers were used for the turbulence 
intensity measurements. Three anemometers were required for the three single sensor probes used. Each 
anemometer was equipped with its own signal conditioner for low-pass filtering, DC offsetting, and 
amplifying. A commercial off the shelf (COTS) 3 -channel 12-bit analog to digital converter with an input 
range of ±5 V was used to digitize the signal conditioned anemometer bridge voltages. A COTS personal 
computer with COTS software was used to control the data acquisition process and to store the digitized 
data on computer hard disk. 



Figure 13. — Layout of the RTD probes mounted on the D corner turning vanes. 
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The hot wire probes were calibrated in place using the velocities measured by the flow angle pressure 
probes. Data was acquired at small enough velocity increments to generate calibration curves for the hot 
wires. The probes were generally used from 50 mph to 175 mph, i.e., within the incompressible range. 
Over this speed range, the bridge output voltage varied between 1 .3 and 1.8 V. 

Test Procedures and Execution 

This section and the subsequent subsections describe aspects of the test planning, setup and execution 
that were critical to the success of the aero-thermal calibration in the IRT. As indicated in Reference 2, it 
is important that tunnel calibration tests are given proper attention to detail and with respect to project 
management. Also, customer requirements must be considered when developing the tunnel calibration 
test plan, test matrix and hardware. This specific test considered needs from both internal and external 
customers requiring improved total temperature measurement accuracy, correspondingly improving total 
temperature calibration. These requirements had a direct impact on the development of the new test 
hardware as well as the test planning and execution. 

Test Matrix 

The recommendations made following the 2000 aero-thermal calibration (Ref. 4) were implemented 
for the 2004 full aero-thermal calibration. The “aero” and “thermal” data collection were separated; the 
9-ft survey rake was used to collect the “aero” data at ambient temperature conditions and the 2-D RTD 
array was used to collect all temperature calibration data. The “aero” test matrix, seen in Table 1 was 
completed in 3 test days. Only 1 day was needed to complete the RTD array test matrix, seen in Table 2. 

Table 2 also shows the update in how the temperature calibration test matrix was constructed. In 
previous tests, the total temperature calibration test matrix was set up to cover the operating range of the 
IRT by selecting test points at minimum and maximum conditions and at regular intervals in between 
based on the D-corner total temperature. Starting in 2004, the temperature calibration test matrix was 
based on the test section static temperature as outlined in the guidelines provided in Reference 1 . The 
actual test points are still the test section airspeed and D-corner total temperature, but the total 
temperature setting is determined from the airspeed and desired test section static temperature. The order 
of test points were arranged from highest to lowest in order to more efficiently use the available test time. 


TABLE L— TEST MATRIX FOR THE IRT AERO-THERMAL 
CALIBRATION USING THE 9-ft SURVEY RAKE (2004 TEST) 


Rake height, 
Y , 
in. 

Total temperature, 

r ^T,da\gi 

°c 

Test section airspeed, 

u ts , 

knots 

Spraybar air 
pressure, 
p . 

1 am 

psig 

36 

Ambient 
(nominally 5 °C) 

0, 43, 65, 87, 100, 130, 152, 173, 200, 217, 240, 260, 282, 304, 339 

0, 70 

6, 12, 18, 24, 30, 
42, 48, 54, 60, 66 

ambient 

0, 43, 87, 100, 130, 173, 200, 217, 260, 304 

0 
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TABLE 2.— TEMPERATURE CALIBRATION TEST MATRIX USING THE 2-D RTD ARRAY (2004 TEST). 


Test section 
static 

temperature, 

A, to? 

°c 

Test 

section 

airspeed, 

u ts , 

knots 

D-comer total 
temperature, 

r ^T,davg'> 

°c 

Spraybar 

air 

pressure, 

p . 

1 am 

psig 

Test section 
static 

temperature, 

A, to? 

°c 

Test 

section 

airspeed, 

u ts , 

knots 

D-comer total 
temperature, 

^T,davgi 

°c 

Spraybar 

air 

pressure, 

p . 

1 am 

psig 

-30 

43 

-29.8 

0, 70 

-5.6 

43 

-5.3 

0, 70 

-30 

100 

-28.7 

0, 70 

-5.6 

100 

-4.2 

0, 70 

-30 

130 

-27.8 

0, 70 

-5.6 

130 

-3.3 

0, 70 

-30 

200 

-24.5 

0, 70 

-5.6 

200 

-0.3 

0, 70 

-30 

217 

-23.8 

0, 70 

-5.6 

217 

0.7 

0, 70 

-30 

304 

-17.8 

0, 70 

0 

43 

0.2 

0, 70 

-17.8 

43 

-17.6 

0, 70 

0 

100 

1.3 

0, 70 

-17.8 

100 

-16.4 

0, 70 

0 

130 

2.2 

0, 70 

-17.8 

130 

-15.6 

0, 70 

-0.6 

200 

4.7 

0, 70 

-17.8 

200 

-12.5 

0, 70 

-1.7 

217 

4.6 

0, 70 

-17.8 

217 

-11.6 

0, 70 

-7.2 

304 

4.9 

0, 70 

-17.8 

304 

-5.6 

0, 70 

4.4 

43 

4.7 

0, 70 


3.3 

100 

4.7 

0, 70 


One area where the test matrices deviate from the recommendations in Reference 1 was at the higher 
airspeed and temperature settings because it was not possible to set static temperatures of 4.4 °C at 174 
and 260 knots. The total temperatures were 8.3 and 13.3 °C, respectively, since temperature control is not 
available at total temperature settings above about 4 °C). Therefore, the total temperature (D-corner 
temperature) was set at 5 °C with the lower related static temperatures in the test section as listed in 
Table 2. 

Test points were added to the 2-D RTD Array test matrix to cover static temperature conditions 
around the freezing point, 0 °C, since this is an important operating range for the icing community. Test 
points were also added to the test matrices corresponding to airspeed settings of 100, 200 and 300 knots 
so that flow quality information would be available at these specific settings. 

Test Procedures 

For the 2004 test entry, testing using the 9-ft survey rake was completed in May. Testing using the 
2-D RTD array was delayed until July due to additional time allowed for icing cloud calibration. 

The first test configuration was with the 9-ft survey rake positioned at the test section centerline. Prior 
to the first tunnel run, the flow angle probe, pitch and yaw offset angles were measured with respect to the 
tunnel centerline using a six-degree of freedom digitizing arm. Providing a check, in addition to these 
measurements pitch and yaw offset angles were verified with a digital level, tape measure or ruler. The 
combination of these measurements formed the baseline pitch and yaw offset angles. For each change in 
Z, vertical rake height change, the flow angle probe pitch and yaw offset angles were measured again 
with the digital level and tape measure and compared to the baseline values. All changes from the 
baseline were recorded and were later used to correct the pitch and yaw angles generated by the data 
reduction process. Table 3 shows the typical test conditions encountered during the 9-ft survey rake 
testing. The values listed in Table 3 are the average of all readings collected for each nominal setting 
across all test runs. It is important to note that all testing using the 9-ft survey rake was conducted at 
temperatures of approximately 4 °C {T T j vn > 4 °C). This is done to avoid having the tips of the pressure 
probes freeze during tunnel operation and to decrease the possibility of getting water in the 
instrumentation lines. 
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TABLE 3.— TYPICAL TEST CONDITIONS IN THE IRT DURING THE 9-ft SURVEY 
RAKE TESTING AS PART OF THE 2004 AERO-THERMAL CALIBRATION 


Test 

section 

airspeed, 

u ts 

(knots) 

Spray bar 
air 

pressure, 
p . 

1 air 

(psig) 

Bellmouth 

Mach 

number, 

Mbm 

Bellmouth 

total 

pressure, 

E T,bm 

(psia) 

Bellmouth 

static 

pressure, 

Ps,bm 

(psia) 

D-comer average total 
temperature, 

r ^T,davg 

(°R) (°C) 

Drive 
system 
fan speed 
(rpm) 

43.5 

0 

0.067 

14.308 

14.264 

510.3 

5.4 

76 

51.3 

70 

0.079 

14.315 

14.252 

502.7 

6.1 

76 

64.5 

0 

0.099 

14.310 

14.213 

505.7 

7.8 

109 

70.1 

70 

0.107 

14.314 

14.200 

507.1 

8.6 

109 

87.3 

0 

0.134 

14.308 

14.131 

507.9 

9.0 

143 

90.4 

70 

0.139 

14.311 

14.120 

504.1 

6.9 

143 

99.4 

0 

0.153 

14.305 

14.074 

504.7 

7.2 

163 

101.6 

70 

0.156 

14.310 

14.070 

505.7 

7.8 

163 

131.0 

0 

0.201 

14.299 

13.902 

507.5 

8.8 

212 

131.5 

70 

0.202 

14.304 

13.905 

508.0 

9.0 

212 

150.9 

0 

0.232 

14.296 

13.770 

505.8 

7.9 

242 

150.6 

5 

0.232 

14.300 

13.775 

506.0 

8.0 

242 

151.1 

20 

0.233 

14.301 

13.771 

503.9 

6.8 

242 

152.1 

70 

0.234 

14.301 

13.768 

507.3 

8.7 

242 

172.2 

0 

0.265 

14.299 

13.618 

507.2 

8.6 

275 

173.4 

70 

0.268 

14.306 

13.611 

503.9 

6.8 

275 

198.8 

0 

0.306 

14.290 

13.394 

510.5 

10.4 

315 

200.1 

70 

0.309 

14.295 

13.381 

507.5 

8.8 

315 

216.4 

0 

0.334 

14.290 

13.230 

509.1 

9.7 

342 

218.6 

70 

0.337 

14.293 

13.212 

509.3 

9.8 

342 

237.3 

0 

0.365 

14.291 

13.037 

515.1 

13.0 

373 

237.9 

70 

0.366 

14.294 

13.032 

514.4 

12.6 

373 

257.6 

0 

0.397 

14.284 

12.817 

515.4 

13.2 

404 

259.0 

70 

0.399 

14.287 

12.803 

515.0 

13.0 

404 

280.3 

0 

0.430 

14.280 

12.573 

521.1 

16.4 

439 

281.8 

70 

0.434 

14.283 

12.552 

519.0 

15.2 

439 

300.9 

0 

0.466 

14.270 

12.296 

515.2 

13.1 

473 

302.9 

70 

0.468 

14.271 

12.283 

518.0 

14.6 

472 


The Escort system was used to collect data for both the 9-ft survey rake and the 2-D RTD array. For 
each test condition, two to three readings were taken at a scan rate of one per second, for 20 sec, and 
averaged to generate the final reading. 

All three probe types were recorded simultaneously by individual analog-to-digital converter data 
acquisition boards integrated into a personal computer. The data acquisition boards were configured to 
sample three or six channels simultaneously, depending on whether a single sensor or dual sensor X probe 
was used. The sensors were set to acquire data at a rate of 10,000 samples per second and 100,000 data 
points per channel. This is the equivalence of 10 sec of data for each channel. The signal conditioners for 
each hot wire probe were configured to a low pass filter at 4000 Hz, DC offset of-1 V, and a gain of 5. 

The addition of the 2-D RTD array increased efficiency of gathering total temperature calibration and 
flow quality data. The array is mounted at the test section test plane and requires minimal documentation 
of axial station and probe measurement prior to operation. Unlike the 9 ft survey rake no detailed 
alignment procedure of the RTD probes is required when using the 2-D RTD array. Checks of the RTD 
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array instrumentation connections and structure are completed prior to the tunnel run. The most efficient 
way to complete the RTD array test matrix is to start at the coldest bellmouth total temperature condition 
and work through the matrix toward the higher conditions. This is done because it is faster and easier to 
change airspeed settings at a given total temperature, rather than hold the airspeed constant and change 
temperature. The order of test points for the RTD array test matrix, seen in Table 2, is listed by increasing 
total temperature. 


Data Reduction and Analysis 

This section provides a detailed description of the steps taken in the reduction and analysis of all data 
collected using the 9-ft survey rake and the 2-D RTD array. Included in this section are instrumentation 
calibrations, data corrections and conversion equations. The conversion equations were used to convert 
the measured pressures and temperatures into Mach number, flow angle and turbulence data. For all of 
the equations listed, the specific heat ratio, y, is 1 .4 and the specific gas constant for air, R , is equal to 
1716 ft 2 /(sec 2 °R). 

U.S. customary units of measurement are used in the following data reduction. However, the icing 
community prefers a mixed set of units including airspeed in knots and temperature in degrees Celsius 
(°C). Therefore, the details on the data reduction are presented as they are carried out in English units. 
The final results have been converted to and are presented in the units used by the icing community. 

Facility Instrumentation 

The measured and calculated facility parameters that are of greatest importance to the IRT aero- 
thermal calibration are the bellmouth total pressure, bellmouth static pressure, and the D-corner total 
temperature. The bellmouth total pressure is the average of the two measured total pressures from the 
north wall and south wall 


r T,bm 


= {Pi 


T , north 


+ p, 


T, south 


)/2 


In the IRT, the total pressure and differential pressures are measured. The static pressure for each 
bellmouth probe is calculated from their difference. The north and south probe are each calculated using 
the equation 


p s =p t -ap 

The bellmouth static pressure is again the average of the two probe calculations 

Ps,bm — (^S , north ^S, south )/^ 

The facility total temperature, T T ^ avg is the average of the 24 RTDs mounted on the leading edge of the D- 
corner turning vanes 


1 24 
-Yt t 

7 A ^ 1 


L T , davg ^ , / PTJJ 
i=l 


In addition to the basic parameters, some calculated parameters were also needed based on the bellmouth 
total and static pressure measurements, Pr,bm and P s ,bm , respectively. The calculated bellmouth Mach 
number is 
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M bm 


Y-l 


f p A 

r T,bm 

\^Ps,bm J 


y-l 


-1 


The dynamic pressure based on bellmouth conditions is 


Qbm ^ ' ^S,bm ' ^ bi 


bm 


9-ft Survey Rake Calculations 

Data collected by the 9-ft survey rake is used to construct calibration curves for total pressure and 
Mach number. Additionally, the collected data provides flow quality information on total pressure, static 
pressure, Mach number, flow angularity and turbulence. For the following calculations, subscript “bm” 
refers to measured or calculated parameters associated with the two bellmouth pitot-static probes located 
on the inner and outer tunnel walls at the exit of the bellmouth, as shown in Figure 4. The subscript 
66 rake ” refers to measured or calculated parameters associated with any of the probes on the 9-ft survey 
rake. Additionally, the subscript 66 local” refers to the measured or calculated corrected “rake” parameters, 
which have been corrected by individual probe calibration coefficients and represent the true “local” 
properties in the test section. Note that all “local” total pressures, static pressures, and Mach numbers are 
normalized by bellmouth (“bm”) parameters to arrive at recovery ratios. All of the compressible flow 
equations used in the data reduction can be found in Reference 9. 

Basic Flow Parameters 

The average static pressure sensed by each pressure probe was calculated by averaging the four static 
pressure taps on the probe body 


p (^6 + Pi + ^8 + P 9 ) 

T S,rake ^ 

The probe total pressure is simply the pressure measured by port P 5 (Pr,rake = Ps)- Mach number is 
calculated using the ratio of the static to total pressure measured by each probe 


M rake 


2 

y-l 


y-l 

f n \ 

Pr,rake Y 

p 

\ 1 S,rake J 


-1 


The correction equations for total and static pressure are based on the rake Mach number. C Q and C q are 
functions of M rake and are experimentally determined during calibration of the pressure probe flow angle. 
Values for these coefficients can be found in Reference 8. 


Pt, local Pt , rake rake ) ’ \Pt , rake Ps , rake ] 

( Pf , rake ~ Ps , rake ) 


r S, local r T , local 


Cq rake ) 


The local test section Mach number for each probe is determined using the Pt, heal and Ps, local calculated 
above 
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Local dynamic pressure in the test section is then determined 



The thermocouples on the 9-ft survey rake are only used to monitor the temperature during the aero 
portion of the calibration. Total temperature calibration and flow quality data are collected using the 2-D 
RTD array. The following equations are used to reduce the collected temperature data. The equations 
containing 46 local” and “rake” parameters apply to all eleven individual total temperature probes mounted 
on the 9-ft survey rake. The corrected local total temperature, T T j ocal , is calculated as a function of T T?rake 
and Mi^i 


?T, local Cf local ) * ,rake 


C, is the recovery coefficient for each individual total temperature probe and is a function of M /oca /. These 
coefficients were determined experimentally when the probes were calibrated in the NASA Glenn 3. 5 -in. 
diameter free jet calibration facility. The total temperature probes gradually predict temperatures lower 
than the actual temperature as the local Mach number increases. For example, at a local Mach number of 
0.1, typical values of C t are 1.001 and at a local Mach number of 0.6, the values of C t are generally 
around 1.005. The values of C t are listed in Reference 10. The local static temperature was calculated to 
monitor data and conditions during the test 



The local airspeed is based on a matched set of one pressure probe and one thermocouple from the 9- 
ft survey rake. Each pressure probe was matched to the thermocouple mounted directly under it on the 
rake. The local airspeed for each set of probes was calculated using the following equation 


Flow Angularity 

The equations in this section were used to generate pitch and yaw flow angle data from the raw 
pressure data obtained from the flow angle pressure probes (recall Figure 6). The average of the four flow 
angle pressure measurements is used as part of the calculation for both pitch and yaw flow angles 


p = 

1 avg 


P\ + Py + P\ + Pa 


Pressure coefficients are calculated for both the pitch angle and yaw angle based on the five pressure 
ports on the hemispherical head of each probe. The pitch angle is related to the pressure difference 
between the ports in the pitch plane, ports P\ and P 3 , while the yaw angle is related to the difference 
between the yaw plane ports, P 2 and P 4 , as shown in the following equations 
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c 


( P 5 ~ Evg ) 


and Cp = 


{Pa-P 0 

(^5 ~ P avg) 


The pitch angle, a, and yaw angle, P, are calculated as a function of M rake , C a and Cp. 

^ — 0,a rake) A ^l,a rake) ' C a ,arm ^ a , lev el 

P — ^0, P m ke) A ^2, (3 rake) ’ ^(3 ^(3, arm ^f3,rw/er 

Pressures Pi through P 5 are measured directly from the flow angle pressure probes. The coefficients 
^o,a ? ^i,a , Ko$ , and K 2 , p were experimentally determined when the flow angle pressure probes were 
calibrated in the NASA Glenn 3. 5 -in. diameter free jet calibration facility and are all functions of M ra ke- 
Values for these coefficients are documented in Reference 8. The A a arm and Ap tarm values were determined 
by measuring the pitch and yaw offset angles for the probes with a 6 degree of freedom digitizing arm and 
were constant for every Z elevation. The A a j eV ei and Ap ?rw/er values were equal to 0° for the Z=36-in. 
elevation since this elevation is where the baseline measurements were taken and were different for every 
other Z elevation. These values were determined by taking measurements with a digital level and a ruler 
and comparing these measurements to the baseline values. Deviations from the baseline numbers were 
entered as A a, level and Ap ?rM / er . 

Turbulence Intensity 

As was previously indicated, the hot wire probes were calibrated in situ in the IRT. This is because it 
is difficult to match the temperature and static pressure conditions of the IRT with those in the probe 
calibration facility. Data were acquired so that calibration curves could be constructed after acquisition 
was completed. Time averaged bridge voltages were matched to time averaged airspeeds from the flow 
angle pressure probes from the 9-ft survey rake for the same condition. Calibration curves, exponential 
functions with effective air velocity as a function of hot wire bridge voltage, were generated for each hot 
wire probe sensor. The equations used in this curve-fit process are from Reference 1 1 3 and are given 
below, 


U HW,eff,i - U local * C0S ®HW,i i tt _ a r , A 0 A lt -E HW t 

anQ U HW,eff,i ~ A 0,i & HW,i + A \,i e ’ ’ 

where i = 1 for a single, straight-wire sensor and i = 1 to 2 for a dual cross-wire sensor (X probe). Ohwj is 
the angle of the hot wire sensor relative to the freestream velocity vector minus 90°. For a single sensor 
probe, the sensor is normal to the velocity vector and Gnwr , i is (90° - 90° =) 0° for this case. For a dual 
sensor X probe at 0° angle of attack, one sensor is at 45° and the other is at 135° relative to the velocity 
vector. For these two sensors, G H wr,i is (45° - 90° =) ^15° and Ghw,2 is (135° - 90° =) +45°. Ui oca i is the 
freestream velocity, and U H w,eff,i is the velocity component normal to the hot-wire sensor. E HW ,i is the hot- 
wire voltage measured by the anemometer. Once the values of U H w,eff,i are computed, these values along 
with the corresponding values of hot wire anemometer bridge voltage, E HW j can be used to arrive at 
values ofZ 0 ,/, Ay, andZ 2 ,/ via non-linear least squares curve-fitting using the exponential equation given 
above. 


3 


The procedure was changed in subsequent calibrations to use the King’s Law equation and single-wire probes. 
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UHW,eff, 1 



U local 


Figure 14. — Schematic depicting terminology on one wire of X probe. 

The reader may be wondering how accurate flow angles and transverse velocities can be computed by 
the dual sensor X probe if no angular calibration data is used. To answer this question, two key 
assumptions are called upon. The first key assumption is that the X probes will only be used to assess the 
unsteadiness in flow angle and will not be used to quantify the absolute flow angle. The resulting flow 
angles and transverse velocities will all be centered around zero and have averages essentially at zero. 

This should work fine in the IRT test section since we know that the flow angles will be small, typically 
±2° or less, and will be centered about zero. The second key assumption is that the X probes used in this 
test behave like other X probes previously calibrated through various flow angles. The results for an X 
probe calibrated through various flow angles show that the calibration map will collapse down to one 
curve when effective velocity is plotted versus hot wire anemometer bridge voltage. This typically holds 
true for angles between ±20°. 

Given the exponential curve fit equation used to relate hot wire effective velocity to hot wire 
anemometer bridge voltage, data reduction can continue with the equations given below for a dual sensor 
hot wire X probe. 


Test section airspeed: U HW = yJuj [W e(f i + Uf [W ^ e(f 2 


Pitch flow angle: a HW - 45° - tan 1 


A 


U HW,eff, 2 
k U HW,eff,l J 


Axial velocity component: u HW = U HW • cos(a# r ) 
Transverse velocity component: v HW = U HW • sin(a// r ) 


For a single sensor hot wire probe, the following equalities are used instead of the four equations 
above. 


U HW ~ U HW - U HW e ff X 

Once u H w is available for a single sensor probe, and u HW and v HW are available for dual sensor X 
probe, statistics necessary for computing turbulence intensity can be calculated using the equations given 
below. The value of N was 100,000. The average axial and transverse air velocity components are 
determined using the following equations. 

y N y N 

u HW,avg = — ~^ u HW,i an( l v HW,avg = ~ HW ,i 
i - 1 i - 1 

The magnitude of the velocity vector, also known as airspeed, is determined by combining the axial and 
transverse velocity components. 
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U HW,avg - 

The standard deviation of the axial and transverse velocity components are calculated using the following 
standard statistical equations. 


Finally, the turbulence intensity for the axial and transverse velocity components are defined as 

TI U = auHW xl00% and TI V = ° VHW xl00% 

U HW,avg U H w,avg 

Since in-situ calibration curves were generated, there was no need to correct for changes in 
temperature or static pressure. The calibration curves took care of these inherent variations. In addition, 
new calibration curves were generated for every velocity sweep. 

2-D RTD Array Calculations 

The total temperature measurements made by the 49 RTD probes on the 2-D RTD array only 
required a correction of flow recovery. The determination of the total temperature flow recovery 
correction was described in a previous section of this report. The resulting correction equation is 

Tt, local = (Q + Ci • M ts + C 2 • M} s )• T t array 

Where T^ array is the total temperature measured by each probe and T T j ocal is the corrected total temperature. 
The correction coefficients are listed in Appendix B. Only 20 of the 49 probes on the array were directly 
calibrated for flow recovery and therefore have specific coefficients. These 20 probes were mounted in 
the interior portion of the array. The remaining 29 RTD probes used the average of the correction 
coefficients as listed in Appendix B. 

Measurement and Data Uncertainty 

Details of the data uncertainty analysis are contained in Appendix A. The uncertainty analysis follows 
the standard procedure set in Reference 12. 

Data Analysis 

Standard test procedure was to collect three data readings for each test condition listed in the test 
matrices for the 9-ft survey rake and the 2-D RTD array. Each data reading was the average of 20 scans 
of data, at a data acquisition rate of 1 scan per second. The data readings were used in two combinations 
for development of calibration curves versus assessment of test section flow quality. 

• For the construction of the calibration curves, the readings were used as they were collected 
(average of 20 scans). So, each test condition would contribute three individual readings to the 
development of the calibration curves. 

• Flow quality distribution data plots were created by averaging the three data readings collected at 
each test condition. For data from the 9-ft survey rake, the data represented at each rake position, 
or height, is the average of three separate data readings collected at one test condition over the 
course of 1 1 tunnel runs (one run for each of the 1 1 vertical rake positions). For the 2-D RTD 
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array, the data is again plotted by height of the RTD probes above the test section floor and each 
data set of the average of three readings collected at a given test setting. However, unlike the 9-ft 
rake, for the array, only one tunnel run was needed to collect all the data for a given test setting. 

• To combine data collected using the 9-ft survey rake over several test runs to show the composite 
flow quality in the test section, the pressure, Mach number and airspeed data were normalized 
using data collected simultaneously by the facility bellmouth rakes. 

• The corrected pitch and yaw flow angle data were converted into flow angle vectors to more 
easily visualize the flow direction. The length of the vector represents the resultant flow angle 
magnitude and the orientation of the vector shows the resultant flow direction. The magnitude 
and orientation were calculated as shown below. 

Vector magnitude, V m V m = ^ a 2 +P 2 


f ()>, if a > 0 and P > 0 


Vector orientation, = tan 1 



and = 1 


1 80 - 4>, if oc>0andp>0 
4> H- 1 80, if oc<0andp<0 


[360 - 4>, if oc<0andp>0 


Discussion of Results 

There are three sets of results described in this section: (1) the aero-thermal calibration curves for 
determining the operating conditions in the IRT, (2) detailed aero-thermal flow quality mapping data for 
the IRT test section and (3) uncertainty estimates for the measured and calculated data. 

Aero-Thermal Calibration Curves 

As indicated, one of the primary goals of this test entry was to update the aero-thermal calibration 
curves used to set operating conditions in the IRT test section. A computer program subroutine was 
created following the 2000 calibration. This subroutine contained the calibration relationships and 
additional calculations to fully describe the test section aero-thermal operating point. This subroutine is 
the IRT Aero-Thermal, and abbreviated as IRT AT. The data collected during the 2004 calibration test 
were used to update the IRT AT. 

Updates Based on 2004 Calibration Test 

The existing calibration curves, based on 2000 calibration data, were compared to calibration curves 
developed from the 2004 data. It was not expected that the 2000 and 2004 curves would match due to the 
changes made to the facility, but large variations were also not anticipated. Figure 15 shows the four 
calibration relationships that were created based on the 2004 data, total pressure, static pressure, Mach 
number and total temperature. The primary calibration relationships were changed after review of the 
2004 calibration curves and comparison to the existing year 2000 calibration curves. This change was 
completed by replacing the total pressure curve with the Mach number calibration curve, due to a reduced 
sensitivity on spraybar air injection effects. The 2000 calibration curves do not include spraybar air 
injection. Spraybar air injection and high temperature test points, above 37 °F or 3 °C, were not included 
in the total temperature calibration curve. The computing requirements for IRTAT were updated per the 
new calibration curves (Revision 4, dated July 29, 2004). 4 The calculations used in IRTAT (Revision 4) 
are described below. 


Additional corrections were made to IRTAT so the final version based on the 2004 data was Revision 5, dated 
August 13, 2004. 
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Figure 15. — Aero-thermal calibration curves based on the 2004 full calibration test 


Inputs to IRTAT 

The measured inputs for determining the calibrated test section aero-thermal conditions are acquired 
from two areas. The first area is the north and south bellmouth, or test section inlet, pitot-static probes, 
and the second area is the 24 D-corner RTD probes. Details and definitions on these inputs and their use 
for determining the calibrated parameters are listed below, with the appropriate units listed in parenthesis. 

Bellmouth total pressure (psia), Ptm-> i s the average of the two measured pressures, one each from the 
north and south bellmouth rakes 


r T,bm 


^ i^T,bm,n ^ T,bm,s ) 
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Likewise, the bellmouth delta-pressure (psid), AP bm , is the average of the two delta-pressure 
transducers 


&Pbm ~ ~ Pbm,n + ^bm,s ) 

The bellmouth static pressure (psia), P s ,bm, is calculated from the measured total and differential pressures. 
First the static pressure for each bellmouth rake is determined 

Ps,bm,n ~ ,bm,n ~ ^^bm,n <^nd Ps,bm,s ~ ,bm,s ~ ^Pbm,s 

The bellmouth static pressure is then the average of the static pressure calculated from each bellmouth 
probe 



to 


S,bm,n 


+ Ps 


S.bm.s 


) 


The facility total temperature is the average of the 24 D-comer total temperature RTD measurements 


(°C), T T ,davg 


24 

T T ,davg = 

i - 1 

The Mach number based on bellmouth conditions, M bm , is calculated from the ratio of the bellmouth 
static to bellmouth total pressure using the standard isentropic flow equation (Eq. 44 from Ref 9 with 

y=i-4) 


M 


bm 


( p v 0 - 2857 

, S,bm 

y ,bm J 


Test Section Conditions; Calibrated Values 

The test section calibrated values, Mach number, static pressure and total temperature, are determined 
from calibration curves based on bellmouth pitot-static probes and D-corner RTD measurements. Figure 
15 shows the above mentioned calibration curves. The relationships for the test section conditions are 
listed in the subsequent section. The calibration relationships were based on data collected using the 9-ft 
survey rake and the RTD Array. The 9-ft survey rake measured total and static pressure which is used to 
compute Mach number and the RTD Array measured the total temperature. The calibrated test section 
values are based on conditions measured at the test section turntable centerline and represent average flow 
field values at the core of the test section (the center 25 measurement locations from both the 
aerodynamic and thermal surveys). The pressure and Mach number calibration curves include data for all 
airspeed and spraybar air pressure settings, however, the total temperature calibration curve based on the 
2004 calibration excludes data at the following conditions: T T ^av g > 3 °C, airspeed< 50 knots, P au > 0 psig. 
The calibrated test section static pressure (psia), P SftS9 is determined from the static pressure calibration 
curve, seen in 

Figure 15(b), based on the least squares curve fit to the data 

P S ts = -0.3770668 1 + 1 .02636653 • P SM 
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The calibrated test section Mach number, M ts , is likewise calculated from the Mach number calibration 
curve, shown in 

Figure 15(c), and best fit line 


M ts = 0.00082798+1.01631617- Af ftjn 

Test section total temperature (°C), T TJs , is determined from the calibration curve shown in 
Figure 15(d) and calculated from the following equation 

T Tts = 0.24698 1 + 0.998554 • T Tdavg 

It is important to note that the units for the calibrated test section conditions are not in a consistent set 
of engineering units. For subsequent calculations, the total temperature was converted from metric to 
English 


T Us (°F) = f T Tjs (°C) ~] + 32 and T rjs (°R) = T rjs (°F) + 459.67 
v 5 


Test Section Conditions; Calculated Values 

The remaining test section aero-thermal operating conditions are calculated based on the values of the 
calibrated parameters (Ps,ts, M ts and T^ ts ). The calculated parameters include the total pressure, static 
temperature, airspeed and Reynolds number. The details of these calculations are included in Appendix C. 

Updates Based on 2005 Interim Calibration 

As described in Appendix D of the IRT and Reference 3, the total temperature calibration was 
improved in 2005. The improved total temperature calibration was incorporated into IRT AT as 
Revision 6 (dated April 11, 2005). These updates are documented in Appendix D. 


Aero-Thermal Flow Quality Data 

In addition to creating the aero-thermal calibration relationships, the data collected using the 9-ft 
survey rake and the 2-D RTD array are also used to construct detailed descriptions of the test section flow 
field. The flow quality data for total and static pressure, total temperature, airspeed and Mach number are 
presented in x-y plots (rather than contour plots) to provide as much detail as possible about the flow 
field. Pitch and yaw flow angle data were combined to create flow vectors that provide a combined flow 
angle magnitude and direction. Note that for most of the data presented, the data collected with the 9-ft 
survey rake at 66-in. above the test section floor (maximum height) was omitted as this data was affected 
by boundary layer (either the rake was in the boundary layer or tripped the boundary layer). Also, in some 
cases, data near the test section walls has been truncated due to the choice of y-axis scale; this was done 
in order to show the details of the flow field in the test section core. The total pressure measurement on 
probe number two (18 in. from the inside test section wall) was inoperative during testing (plugged tube). 
These data were omitted from the flow quality plots and account for the gap in the total pressure and 
Mach number distribution data at the position 18 in. from the wall. As noted earlier, 9-ft survey rake data 
were collected at bellmouth total temperature settings greater than 4 °C (about 40 °F). 

Aero-thermal flow quality goals for an icing tunnel are defined in Reference 1 and are summarized 
here in Table 4. In general, the IRT flow quality is very good in terms of the pressure, temperature and 
Mach number (airspeed) distributions and meets or exceeds the goals listed in Table 4 (as will be 
discussed in subsequent sections). While the IRT does not have any conventional flow field manipulators 
(i.e. turbulence reduction screens or flow straightening honeycomb), the heat exchanger does act a flow 
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manipulator and provides a significant pressure drop that evens the flow field in the tunnel. The flow 
angularity in the IRT test section is rather high compared with an aerodynamic wind tunnel of similar 
operating range. This is due to the lack of flow straightening and also the presence of the spraybars in the 
settling chamber upstream of the test section. The IRT is a classic example of the trade-offs made in 
ground testing; in order to create icing conditions, test section aerodynamic flow quality is somewhat 
sacrificed by the omission of normal flow field manipulators since those devices would collect ice and 
effect operation and icing performance of the tunnel. 


TABLE 4.— AERO-THERMAL FLOW QUALITY GOALS FOR AN ICING TUNNEL TEST SECTION 


Parameter 

Measurement 

Uncertainty 

Test section spatial 
uniformity 

Tunnel centerline 
temporal stability 

Airspeed 

+1% 

±2% 

±2% 

Static air temperature between -30 and +5 °C 

±2 °C 

±1 °C 

±0.5 °C 

Flow angularity 

±0.25° 

±2° 

N/A 

Turbulence (P air = l 0 psig) 

±0.25% 

<2% 

±2% 

Turbulence (P air = 70 psig) 

±0.25% 

<2% 

±2% 


Total Pressure 

Figure 16 shows the total pressure flow quality data over the operating range of the IRT. As described 
above, the total pressure data from the 9-ft survey rake were divided by the corresponding bellmouth total 
pressure to remove the run-to-run variations. The data collected from various rake heights for like 
operating conditions were then combined to provide a representation of the test section flow quality in a 
plane located axially at the centerline of the test section turntable. 

With the spraybar air off (P air = 0 psig), the total pressure data show a very uniform distribution 
across the test section with very little variation in the core of the test section. At the lower airspeed 
settings (£4<150 mph or 130 knots) the total pressure variation throughout the entire surveyed area of the 
test section is 0.003 psia or less. At about 200 mph (174 knots) a very slight gradient trend begins to form 
in the test section, with higher total pressures recorded near the inside test section wall. Worst case is at 
the 350 mph (304 knots) setting, where the total pressure variation is about 0.006 psia. 

Adding spraybar air into the tunnel airflow has a measurable effect on the total pressure distribution 
over the entire airspeed range of the facility. In order to gage the worst case effect, data were collected 
with the spraybar air at 70 psig, although there was some data collected at lower air pressure settings. The 
effects of the spraybar air injection are larger at lower test section airspeeds. The largest effects are also 
seen near the test section boundaries, so the spraybar air injection appears to have an effect on the test 
section boundary layer. The flow in the core of the test section is affected minimally. For example, at 100 
mph (87 knots) the total pressure variation in the central 5-ft core area is about 0.003 psia for both the 0 
and 70 psig spraybar air injection cases. The difference is seen in the 2-ft area near each test section wall 
where larger variations were measured (about 0.010 psia near the outside tunnel wall). Similar results are 
seen at higher test section airspeeds. The spraybar air injection also increases the total pressure gradient 
recorded in the test section at airspeeds above 200 mph. At the 300 mph (261 knots), the gradient 
increases from about 0.003 psia at 0 psig spraybar air to about 0.006 psia at 70 psig spraybar air. 

Overall, the IRT test section flow quality is excellent for the total pressure distributions. The 
magnitude of the variations discussed above is quite small and should not have any significant impact of 
the testing conducted in the IRT. 
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Figure 16. — Continued, (e) U ts = 115 mph (100 knots), P a/r =0 psig. (f) U ts = 115 mph (100 knots), P a ,F70 psig. (g) L/f S =150 mph (130 knots), P a/r = 0 
psig. (h) l/f S =150 mph (130 knots), P a/ r= 70 psig. 
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Figure 16. — Continued, (i) U ts = 175 mph (152 knots), P air = 0 psig. (j) U ts = 175 mph (152 knots), P a/r =5 psig. (k) U ts = 175 mph (152 knots), P a/r = 20 
psig. (I) U t s= 175 mph (152 knots), P a/ >=7 0 psig. 
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Figure 16. — Continued, (q) U ts = 225 mph (195 knots), P air = 0 psig. (r) U ts =225 mph (195 knots), P air = 20 psig (no z=36-in. vertical position data), 
(s) U ts =225 mph (195 knots), P a />=7 0 psig. 
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Figure 16. — Concluded, (t) U ts = 300 mph (261 knots), P a /r= 0 psig. (u) U ts =300 mph (261 knots), P a/ p7 0 psig. (v) U ts =350 mph (304 knots), P a/r = 0 psig. 
(w) U ts =350 mph (304 knots), P a/ >=70 psig. 


Mach Number 

The Mach number ratio distribution data are shown in Figure 17. As was the case with the total 
pressure data, the Mach number measurements from the 9-ft survey rake mounted in the test section were 
normalized by the bellmouth Mach number so that the data could be combined to show the overall flow 
quality for a given test section condition. For each set of Mach number ratio data plots, an approximate 
delta Mach number (A M) was added to provide a gage of the absolute variation in the Mach number 
within the survey plane. Note that the y-axis scale is expanded to show flow field detail. 

The overall Mach number variation for conditions with no spraybar air injection (P air = 0 psig) was 
0.002 to 0.004, neglecting points near the test section surfaces that may have been affected by the 
boundary layer. The shape of the lateral distribution were consistent at each rake vertical position with 
slightly higher Mach numbers recorded nearer the outside test section wall and the variation was greater 
at the higher Mach number settings (note that the difference in Mach number from inside wall to outside 
wall was only 0.001 to 0.003). 

Spraybar air injection did not have a significant effect on the Mach number distributions, particularly 
in the core of the test section. At the lowest airspeed settings, the spraybar air injection appeared to 
increase the boundary layer thickness near the test section walls (as was also seen in the total pressure 
data) resulting in lower Mach number readings at some rake positions for probes near the walls. For 
nominal airspeed settings of 200 mph (174 knots) and higher however, 70 psig spraybar air injection had 
little or no effect on the Mach number distribution data. 

The Mach number distributions at the test section survey plane are very good over the entire 
operating range of the IRT, including those conditions with 70 psig spraybar air injection. Considering 
that there are no flow quality manipulators (i.e. screens or honeycomb) installed in the settling chamber 
upstream of the test section, the Mach number flow quality is remarkably good (the flat heat exchangers 
and D-comer turning vanes do have a positive effect on the test section flow quality, as described in 
Ref. 4). 

Total Temperature 

For the icing tests that are conducted in the IRT, temperature is the most important of the aero- 
thermal measurements. Both the accuracy of the total temperature measurement and the uniformity of 
temperature in the test section are critical. Because of the critical nature of the temperature measurements, 
the IRT staff is always looking at ways to improve the quality of the total temperature data collected in 
the IRT and the overall flow quality in terms of the total temperature distributions in the test section. As 
mentioned above and in Reference 3, both the accuracy of the total temperature measurements and the 
test hardware used to survey the IRT test section were improved prior to the 2004 aero-thermal 
calibration. These improvements include upgrading the facility and calibration total temperature 
instrumentation from thermocouples to RTDs and using a dedicated total temperature survey tool, the 2- 
D RTD array, for test section measurements. 

The total temperature setting in the IRT is controlled by the two flat heat exchangers located in the C- 
D leg of the tunnel loop (Figure 1). The 24 RTDs mounted on the leading edge of the D-comer turning 
vanes (Figure 14) are used to monitor the temperature uniformity of the air exiting the heat exchangers. 
The control system is used to vary the amount of coolant to each of the 8 heat exchanger zones in order to 
produce a uniform total temperature profile in D-corner. The goal is to control the temperature 
distribution to within ±0.5 °C within the D-comer survey plane (this exceeds the goal set forth in 
Ref. 1 and listed in Table 4). Temperature conditions can be accurately set and controlled from about 4 to 
-30 °C. The critical operating point for the IRT is U ts = 230 mph (200 knots) and T T ^ avg = 0 °C. Note that 
there were some points called out in the test matrix (Table 4) that required total temperature settings 
greater than 4 °C to achieve the desired static temperature. These conditions are beyond the effective 
controllable range of the IRT. 
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Figure 17. — Continued, (m) U ts =200 mph, M ts = 0.265, P a/ >=0 psig. (n) U ts = 200 mph, M ts =0.265, P a/r =70 psig. (o) L/ te = 250 mph, M ts =0.332, 
Pair= 0 psig. (p) Ut s =250 mph, M ts =0.332, P a/r =70 psig. 
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Figure 18 shows the resulting uniformity of the D-comer total temperature for the settings used 
during the aero-thermal calibration. Note that for most all conditions below 3 °C the variation is D-comer 
total temperature readings are within or nearly within the goal. Exceptions are at very low temperatures 
and very low speeds where variations greater than ±1 °C where recorded. The critical operating regime 
for the icing community is airspeed from 1 15 to 300 mph (100 to 260 knots) and temperatures from 
freezing to -20 °C; for those settings the D-comer total temperature distribution was consistently within 
the ±0.5 °C goal. 

The test section total temperature data were collected using the 2-D RTD array. The total temperature 
distributions for each test condition are shown in Figure 19. The data in Figure 19 were sorted first by 
airspeed (lowest to highest) then by total temperature (highest to lowest) for each airspeed setting group. 
Each data plot in Figure 19 shows the test section total temperature distribution for a given airspeed, total 
temperature and spraybar air setting; that data are grouped by vertical height within each plot. Note that 
probe number 3 (top row of probes on the array, third from the inside test section wall) was inoperative 
for most of the test. 

The following are general observations of the test section total temperature uniformity: 

• D-comer total temperature settings greater than 3 °C (T T , davg > 3 °C). — These settings are warmer 
than the controlled operating range of the IRT heat exchange and therefore produce a much 
higher variation in the test section and D-corner total temperature distributions as compared to 
settings less than 3 °C (Figure 18 and Figure 19). The D-comer and test section total temperature 
gradient both show higher temperatures near the tunnel ceiling with a gradual decrease in 
temperature toward the floor for all conditions where T T j avg > 3 °C. As seen by the subsequent data 
plots, these “high temperature” settings are not indicative of the IRT total temperature flow 
quality at settings used for icing conditions. 

• D-comer total temperature settings less than 3 °C (T T , davg < 3 °C). — For total temperature settings 
in the effective operating range of the heat exchangers (3 to -30 °C), the total temperature 
variation in the core of the test section is reduced to 1 °C or less (the core of the test section is 
defined as the center 5- by 5-ft area). Many conditions show variation on the order of 0.5 °C 
(peak-to-peak). 

• Nominal test section airspeed of 50 mph (£4=50 mph =43 knots). — At the very low speed 
condition, the spraybar air injection has a significant effect on the test section total temperature 
distribution at all temperature settings. At P a i r = 0 psig, the distributions are uniform and show 
little variation over test section cross-section (horizontal and vertical distributions), generally on 
the order of 1 °C or less. Spraybar air is heated to nominally 80 °C at P a ir = 70 psig there is a 
noticeable horizontal gradient of about 2 °C across the test section (higher temperatures recorded 
near the outside wall) as well as vertical variation of about 2 °C. 

• Temperature spikes near test section boundaries. — There are several instances of temperature 
spikes recorded near the test section floor or ceiling (probe vertical position of 9 and 63 in. above 
the floor, respectively). For example, there is a spike seen near the ceiling, 54 in. from the inside 
wall, at the £4=50 mph (43 knots); worst case is a 2 °C spike at T T ^ davg =-5A °C. Also, there is an 
occasional temperature spike at the floor, generally 72 in. from the inside wall at U ts of 1 15 mph 
(100 knots) and higher). These spikes are near enough to the test section boundaries that they 
should not have any impact on the core flow temperature distributions. 

• Primary operating range (1 15<£4<300 mph, 3<7>^ vg <-20 °C). — For all of the critical operating 
conditions (airspeed and total temperature combinations), the total temperature profiles in the 
core of the test section show very uniform conditions that meet the ±0.5 °C goal. In fact, for many 
conditions, the total temperature variation in the core of the test section is 0.5 °C peak-to-peak for 
both spraybar air {P a ir) off and 70 psig conditions. 

• Critical (freezing) conditions (T T , davg =0 °C). — Temperature conditions at or about the freezing 
point of water are of critical importance in the aircraft icing community for all airspeed settings in 
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the IRT. These total temperature distributions in the IRT test section at freezing conditions are 
excellent; at all airspeed settings the variation in the core of the test section is less than the 
recommended goal of ±1 °C stipulated in Reference 1 and for most cases meets the ±0.5 °C goal 
used in the IRT. 
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Figure 18. — Total temperature distribution along the D-corner turning vanes as measured by the D-corner RTD 
probes. See Figure 13 for location of probes relative to the heat exchanger, (a) T T ,davg>4 °C. (b) -0.5>T T ,davg>-3 °C. 
(c) —3 >Tr,davg > —9 °C. 
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Figure 18. — Concluded, (d) -1 '\>T T ,davg>-22 °C. (e) -22>7>,dai/g>-31 °C. 
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Distance from inside tunnel wall, y, inches Distance from inside tunnel wall, y, inches 

Figure 19. — Continued, (i) U ts =50.0 mph (43.5 knots), T T ,dav g =- 29.7 °C, P a/ >=0 psig. (j) U ts =5 0.2 mph (43.7 knots), T T ,davg=- 29.6 °C, P a/r =7 0 psig. 
(k) L/ te =115.6 mph (100.5 knots), T T ,davg=4.8 °C, P air = 0 psig. (I) L/ te =116.4 mph (101.2 knots), T T ,davg=5- 3 °C, P a/ >=7 0 psig. 
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Distance from inside tunnel wall, y , inches Distance from inside tunnel wall, y , inches 

Figure 19. — Continued, (q) U ts = 116.8 mph (101.5 knots), 7>,dai/g=-16.4 °C, Pair 2 0 psig. (r) U ts = 116.7 mph (101.4 knots), T r ,davg =-16.4 °C, 

Pair 2 70 psig. (s) L/ te =116.6 mph (101.3 knots), 7>,dai/g=-28.4 °C, P a /r=0 psig. (t) L/f S =116.4 mph (101.2 knots), T T ,davg=- 28.6 °C, Pair 2 70 psig. 
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Figure 19. — Continued, (kk) U ts = 254.7 mph (221.3 knots), T T ,davg= 4.9 °C, P ai r= 0 psig. (II) U ts = 254.5 mph (221.2 knots), T T ,davg= 4.8 °C, Pair-7 0 psig. 
(mm) U ts - 253.9 mph (220.6 knots), T T ,davg=0.6 °C, P a irO psig. (nn) U ts =253.7 mph (220.5 knots), T T ,davg=0.7 °C, P a/r =7 0 psig. 
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Flow Angularity 

The test section flow angularity as measured by the 5 -hole probes mounted on the 9-ft survey rake is 
shown in Figure 20. As described above, the 5-hole probes measure the pitch and yaw components of the 
flow angle and are then combined to produce the resultant flow vectors. The vectors shown in Figure 20 
only reflect the magnitude and direction of the local air flow, but not airspeed. For the low speed test 
points (£4=50 mph) the data from rake position 66-in. above the test section floor was deleted due to the 
low speeds measured at that location (the airspeeds were much lower than the calibrated range of the flow 
angle probes). Larger flow angles were consistently recorded near the test section ceiling which may be 
due to the suspected boundary layer interaction noted in Mach number and total pressure distribution data. 

The flow angularity in the IRT test section generally meets the ±2° goal stated in Reference 1 
(Table 4), particularly in the core of the test section (roughly the center 4.5 by 3 ft area). In this central 
core area, the average flow angle vector magnitude is less than 2° for all conditions (including maximum 
spraybar air pressure) as shown in Figure 21. The average flow angle vector magnitude hovers between 
0.7° and 0.9° over the entire airspeed range with no spraybar air injection (74>= 0 psig) This data also 
indicates that there is a significant effect on the average flow angle vector magnitude at maximum 
spraybar air pressure (P air = 70 psig) up to about £4=150 mph. Above £4=150 mph (130 knots), spraybar 
air pressure effect on the average flow angle vector magnitude is less than 0.1 degree between P air = 70 and 
P a ir = 0 psig. The limited data collected at lower spraybar air pressure settings (P flZ >= 5 and 20 psig) show no 
significant difference in average flow angle vector magnitude compared to spraybar air off (P air = 0 psig). 

Overall, the flow angularity in the IRT test section is acceptable for an icing facility and actually is 
not bad for a subsonic wind tunnel that does not have any flow straightening honeycomb in the settling 
chamber. Other than the vortices in the test section comers at low airspeed and high spraybar air settings, 
there are no obvious flow patterns shown in the vector data. Most of the variation in the flow angle 
orientation is probably due to disturbances generated by the spraybars (wakes). 

Turbulence 

Thermal anemometry was used to measure the turbulence levels in the IRT test section. Three hot- 
wire probes were mounted on the 9-ft survey rake. The turbulence data is shown in Figure 22. The IRT 
will not have the low-turbulence levels of a wind tunnel used for aerodynamics testing due to the lack of 
flow manipulators like screens and honeycomb and due to the influence of the spraybars. In general, the 
axial turbulence varies between 0.2 and 1.0-percent with no air injection from the spraybar (P air = 0). The 
traverse turbulence is a little higher, ranging from 0.5 to 2.0-percent. Note also that the turbulence levels 
are higher near the test section centerline due to the wake from the spraybar centerline support. Higher 
turbulence levels were also recorded near the test section floor and ceiling, possible due to boundary layer 
effects. The injection of spraybar air has a significant effect on the test section turbulence, particularly at 
the lower airspeed settings. 
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Figure 20. — Flow angle data from the IRT test section collected using the 9-ft survey rake mounted at 1 1 vertical positions over several tunnel runs. 
Nominal conditions are used to denote each data set. (a) U ts =5 0 mph (43 knots), P a /r=0 psig. (b) U ts =5 0 mph (43 knots), P a />=7 0 psig. (c) L/ te =100 mph 
(87 knots), P ai r= 0 psig. (d) L/ te =100 mph (87 knots), P a/r =7 0 psig. 
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inside wall, probe 2 was 53.7 in. and probe 3 was 76.2 in. (a) Test section airspeed of 50 mph 
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Figure 22. — Continued, (b) 100 mph 
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Figure 22. — Continued, (c) 115 mph 
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Figure 22. — Continued, (d) 150 mph 




* 

o 

in 


* 

in 


|J Qi -^JQ *Jp 

gf* 

o in o m o in 

^ cn d e\i csi ^ 

(%) Aiisuajui aau3|nqjnj[ 


□ 


# 

in 

o 


o 

d 


NAS A/CR— 20 1 0-2 1 6733 


67 


Figure 22. — Continued, (e) 175 mph, 0 and 70 psi. 
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Figure 22. — Continued, (f) 175 mph, 5 and 20 psi. 
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Figure 22. — Concluded, (g) 200 mph. 


Summary of Test Results 

Aero-thermal calibration testing was completed following facility modifications to the NASA Glenn 
Icing Research Tunnel. The original inlet guide vanes upstream of the IRT drive fan were replaced due to 
the age and structural concerns. The new inlet guide vanes were designed to provide more efficient 
operation of the IRT. During January 2004 a full aero-thermal calibration was completed to fully assess 
the impact of the new inlet guide vanes on the IRT test section calibration and flow quality. An interim 
calibration was completed in January 2005. The following were the primary accomplishments and results 
from the 2004 to 2005 aero-thermal calibration test entries in the IRT. 

From the January 2004 full aero-thermal calibration: 

• Completed the planned test matrix that covered the entire operating of the IRT (airspeed, 
temperature and spraybar air injection). 

• From the collected data set, a new set of aero-thermal calibration curves were constructed. These 
calibration curves provided an accurate characterization of the IRT aero-thermal test section 
conditions over the operation range of the facility. 

• The data set also provided a detailed flow field map of the aero-thermal flow quality in the IRT 
test section. 

In addition to providing a general check of the IRT aero-thermal operation and calibration, the data 
collected during the January 2005 interim aero-thermal calibration provided an improved total 
temperature calibration model. 


Recommendations 

The following recommendations were developed following each calibration test entry. Some of the 
recommendations from the 2004 test were implemented during the 2005 test. 

• Following the 2004 test, it was noted that the aero-thermal calibration curves should be refined to 
better model the affects of the spraybar air. This was implemented in the calibration curves 
developed following the 2005 interim calibration. 

• During the 2004 calibration, there some test conditions that required additional stabilization time, 
particularly in terms of total temperature. It was recommended that a means of monitoring the 
stability of the total temperature setting was needed to ensure steady test conditions. The facility 
electrical engineers developed a graphical representation that monitored the real-time, running 
average and standard deviation in the D-corner temperature measurements to monitor data 
stability. This was implemented for the 2005 testing. 

• The data collected in 2004 and 2005 showed that spraybar air injected into the airstream has a 
very noticeable effect on the test section operating conditions as the low speed settings (less than 
100 knots). Additional testing is needed to determine (1) the test section airspeed at which 
spraybar air pressure no longer has a significant effect and (2) the minimum spraybar air pressure 
that causes a noticeable effect at the low speed settings. 
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Appendix A. — Uncertainty Estimates 

This analysis was conducted based on the AIAA Standard “Assessment of Experimental Uncertainty 
with Application to Wind Tunnel Testing” and its companion Guide, References 12 and 13. 

A.l General Uncertainty Analysis 

The accuracy of a measurement is the true value in relation to the measured value. In 
experimentation, the true value is unknown. Therefore, the uncertainty of the measured value is 
determined using the bias and precision errors. The bias error is the fixed or constant error that occurs 
with each measurement, while the precision error is the random or repeatability error. 

General uncertainty analysis (Ref. 14) uses the experimental result, r, with variables up to Xj as seen 
in Equation (Al) 


r = r(X l9 X 29 ...JCj) (Al) 

The general uncertainty equation takes the form seen in Equation (A2) with the expanded form shown as 


-.1/ 


U r = 


2AA 


(A2) 


. 1=1 




u = 


dr 

JX 


u 


x , 


+ 


dr 

aV 




+ ...- 


d Z u 

[dXj x \ 


(A3) 


In this analysis, the uncertainties of the velocity, Mach number, and flow angularity were determined with 
respect to total temperature, total pressure and static pressure. Additionally, the pressure measurement 
uncertainty was determined with respect to the Escort D+ data acquisition system for the IRT facility, 
including the pressure transducer and amplifier. 


A.2 Hardware RTD Temperature Uncertainty 

The individual component uncertainties were reviewed for the RTD total temperature probe system. 
Uncertainty values are given in Table A.l. The component uncertainties were calculated using a root- 
sum-squared (RSS) technique which is a method of determining errors. 

RSS = Ja 2 +b 2 +c 2 (A4) 

The component uncertainties calculated using the RSS technique yielded a system uncertainty of 
±0.10 °C for a probe that had been individually calibrated for flow recovery and ±0.16 °C for a probe that 
had not been individually calibrated for flow recovery. The uncertainties associated with the RTD probes 
are more than an order of magnitude improvement over the uncertainty estimate for the thermocouples on 
the 9-ft survey rake. The system uncertainty associated with the 9-ft survey rake is ±1.1 °C. 4,5 The RSS 
computations for the RTD probes are show below: 

For probes that were tested for total temperature flow recovery: 

RSS = V(0.030 2 + 0.080 2 + 0.006 2 + 0.026 2 + 0.040 2 + 0.030 2 ) = ±0.10 °C 
For probes using average flow recovery coefficients: 
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RSS = V(0.030 2 + 0.080 2 + 0.006 2 + 0.026 2 + 0.040 2 + 0.123 2 ) = ±0.16 °C 


TABLE A.I.— RTD TOTAL TEMPERATURE PROBE SYSTEM COMPONENT UNCERTAINTIES. 


Component 

uncertainty, 

°C 

Component description 

±0.030 

RTD capsule, 0.01% DIN curve 

±0.080 

RTD temperature transmitter 

±0.006 

Resistor to convert transmitter 4 to 20 mA output to DC V (USF 240 250 £2 0.01% 2 ppm/°C) 

±0.026 

Escort data system analog to digital conversion error 

±0.040 

Flow recovery reference temperature measurement error (jet plenum reference) 

±0.030 

Flow recovery curve uncertainty (individual probe cal) 

±0.123 

Flow recovery curve uncertainty (probe not calibrated for flow recovery, average coefficients) 


A.3 Hardware Pressure Uncertainty 
A.3.1 Total Pressure 

Two pitot-static probes, one mounted on the IRT test section vertical north wall and the other on the 
vertical south wall, are used to measure total and differential pressure. The total pressure is measured 
using a 0- to 20-psia Setra Model 270 transducer (Setra Systems, Inc., Boxborough, MA), and the 
differential pressure is measured using a 0- to 3-psid Setra Model 239 transducer. 

The IRT facility data acquisition system utilized the ESCORT D± with 470 Neff amplifier system 
and the Setra Model 270 pressure transducer for the total pressure measurement. 

The calculations shown below detail the bias error of the system. 

B = ±VS1 2 + B2 2 + R3 2 + R4 2 + B5 2 + R6 2 + 57 2 + R8 2 + R9 2 
B = ±7(0.002)2 + (0.010) 2 + (0.001) 2 + (0.001) 2 + (0.010) 2 + (0.004) 2 + (0.001) 2 

B = ±0.0149 psia 


TABLE A.2. — TOTAL PRESSURE COMPONENT BIAS UNCERTAINTIES 


Component uncertainty, 
psia 

Component description 

±0.002 

Cal Lab calibration error 

±0.010 

Non-linearity 

±0.001 

Thermal Zero shift 

±0.001 

Thermal span shift 

±0.010 

Neff gain accuracy 

±0.004 

Neff non-linearity 

±0.000 

Channel-to-channel offset 

±0.000 

Common mode voltage 

±0.001 

Digitizing Error 


Precision error 


5 = +±S1 2 + S2~ 2 + S3 2 + S4 2 + S5 2 + S6 2 
S = ±V( 0.002) 2 + (0.002) 2 + (0.001) 2 + (0.004) 2 + (0.003) 2 + (0.002) 2 

S = ±0.0117 psia 
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TABLE A.3.— ' TOTAL PRESSURE COMPONENT PRECISION UNCERTAINTIES 


Component uncertainty, 
psia 

Component description 

±0.002 

Repeatability 

±0.002 

Hysteresis 

±0.001 

Excitation voltage effect 

±0.004 

Neff Gain Stability 

±0.003 

Zero Stability 

±0.002 

Noise 


Total pressure measurement uncertainty 


Urss 



' 2 S V 


(A5) 


Urss = U Pt =±0.027741 psia 


U % = ( Urss/Pr ) = (U Pt /P T ) = 0.001387 = 0.13% 


A.3.2 Dynamic Pressure 

Equation (A6) was found using the Ames Tables (4) and assuming an isentropically perfect gas. 
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A.3.3 Static Pressure 


Using the Ames Tables (4) and assuming an isentropically perfect gas. 
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(A 12) 


(A13) 


(A 14) 


(A15) 


A.4 Airspeed Calculations 

Using the velocity equation, Equation (A 16), the speed of sound equation, Equation (A 18), and the 
isentropic identity shown, Equation (A 19), the velocity was solved for 


V = M -a 


M = 


f 2 1 


f p 0 

U-u 


UJ 




a - pRT 



(A 16) 
(A 17) 

(A 18) 
(A 19) 
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2 P T 

— y Rt t tt 
Y-U V p s. 


The general expression for uncertainty is given by Equation (A21). U Tt , is the uncertainty in total 
temperature, U P , is the uncertainty in static pressure, and U P , is the uncertainty in total pressure. 


U v = —U Tt + —U Ps + —Up, 
\dT T Tt ldP s Ps \dP T Pt 


To produce algebraic simplification, the uncertainty expression is divided by V 2 and shown below. 


j v f r i ev TT i 2 r 1 8 v TT i 2 r i sv TT 

— = + Up + Up 

V V dTr 1 V dP, s V dPr 1 


Taking the partial derivative of the velocity equation with respect to total temperature, T t , static 
pressure, P s , and total pressure, P T , the uncertainty of velocity is seen in Equation (A26). 
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A.5 Mach Number Calculations 

Using the MuPAD feature in the symbolic toolbox in Matlab the partial derivatives for Mach number 
were solved for and are shown in Equation (A27) and (A28). 
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(A28) 


The partial derivatives were then plugged back into Equation (A3) yielding the following for the 
uncertainty for Mach. 
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(A29) 


A.6 Flow Angle Calculations 

The flow angle calculations are completed utilizing the 9-ft rake flow angularity probes. Figure 6 in 
the main body of the document illustrates the locations and descriptions of each of the pressure orifices. 
Equations (A31) and (A3 2) are the equations used to calculate pitch and yaw utilizing the flow angularity 
probes. 


r l-4 ,avg 




(A30) 
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a = KOA + KIA , ^ ^ v + POFF + DP OFF 

V5 ~ P\-A ,avg ) 

(3 = KOB + K2B ( P 4~ p 2) + YOFF+ DYOFF 

y~5 ~ P\-A,avg ) 


(A31) 


(A32) 


The partials of the pitch and yaw equations were taken with respect to PI through P5 and plugged 
into the uncertainty equations. Equation (A33) is the uncertainty equation for pitch, and Equation (A34) is 
the uncertainty equation for yaw. 
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(A34) 


Utilizing the following inputs from the 2004 data set RDG 1530, the following uncertainties were 
determined. The inputs and the uncertainties calculated for the single data point are shown in Table A.4. 


TABLE A.4.— FACILITY UNCERTAINTIES 


| Input 

Uncertainties | 

Total Pressure (psi) 

14.315 

Total Pressure (psia) 

0.027741 

Static Pressure (psi) 

12.258 

Total Temperature (°C/R) 

0.16/0.08889 

KOA 

-0.1682 

Static Pressure (psia) 

0.0411 

KIA 

12.1366 

Mach Number 

0.0061 

KOB 

0.5052 

Dynamic Pressure 

0.0582 

K2B 

11.9173 

Velocity (ft/sec) 

6.7359 



Pitch (degrees) 

0.0435 



Yaw (degrees) 

0.0427 
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Appendix B. — RTD Probe Flow Recovery Coefficients 

Table B.l contains the total temperature flow recovery curve fit coefficients for each of the RTD 
probes used on the 2-D RTD Array. The corrected total temperature measurement for a given probe is 
determined from the following equation 

Tt, local = (O) + C\ * M ts + C 2 • T t .array 

where T TArray is the total temperature measured by the RTD probe and T T j oca i is the corrected local total 
temperature. 


TABLE B.L — FLOW RECOVERY CORRECTION COEFFICIENTS FOR 


THE RTD PROBES MOUNTED ON THE 2-D RTD ARRA. 


Probe Serial 
Number 

c 2 

Ci 

C 0 

1 

5.681 198xl0“ 3 

8.859325x10" 

1.0 

2 

7.608867xl(T 3 

5.487745x10" 

1.0 

3 

6.518485xl(T 3 

2.3051 13x10" 

1.0 

4 

6.562552x1 0“ 3 

9.399041x10" 

1.0 

5 

5.788899xl0“ 3 

6.373217x10" 

1.0 

6 

6.482487x1 0“ 3 

1.749681x10" 

1.0 

7 

1.011 49 lxl 0 2 

-1.646186x10" 

1.0 

8 

5.661808x10 3 

5.089893x10" 

1.0 

9 

6.506136x10" 

4.493981x10" 

1.0 

10 

5.600892x10" 

5.937824x10" 

1.0 

11 

6.159188xl0“ 3 

7.799326x10" 

1.0 

12 

4.848719x10" 

1.197796x10" 

1.0 

13 

5.634439x10 3 

1.523486x10" 

1.0 

14 

5.412994x10" 

1.068442xl0“ 3 

1.0 

15 

6.099340x1 0" 

6.670401x10" 

1.0 

16 

5.39961 1x10" 

6.660572x10" 

1.0 

17 

6.631578x10" 

6.137670x10" 

1.0 

18 

6.718590x10" 

6.487126x10" 

1.0 

19 

4.910175x10" 

8.683076x10" 

1.0 

20 

6.384046xl0“ 3 

6.479808x10" 

1.0 

21 

6.823 889x10" 

3.565925x10" 

1.0 

22 

5.007236x10 3 

1.218866x10" 

1.0 

23 

6.436245x10" 

5.701839x10" 

1.0 

Average 

6.881316xl0“ 3 

1.490507x10" 

1.0 
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Appendix C. — Calculated Test Section Conditions 

The following set of equations is used to compute the test section flow conditions based on the 
calibrated test section conditions. The calibrated test section conditions are static pressure, P SitS9 Mach 
number, M ts and total temperature, T Tfts - For all calculations, y=1.4 andi?=1716 ft-lb/(slug-°R). English 
units are used for all calculations; however some units conversions are required to provide operating 
conditions in the units desired by the aircraft icing community. 

The first parameter calculated is the test section total pressure, P T p- The ratio of static to total 
pressure is determined from calibrated Mach number found in Reference 9 as Equation 44. 


rat,ts 


r SAs 


T,ts 




- 7/2 


V 


J 


The total pressure is then calculated using the pressure ratio and the calibrated static pressure. Note that 
the units used in the calculation are psf. 


Pr,ts - 


f p ^ 

r T,ts 
V P S,ts J 


P S,ts - 


S,ts 


rat,ts 


The test section static temperature, is calculated, in °R, as shown below using Reference 9, 
Equation 43 


Tsjs ~ ' 


L T,ts 


1 + 0.2 -Ml 


Test section airspeed U ts in ft/sec, is determined from the definition of Mach number, M=V/a. Solving 
for V and substituting the definition for the speed of sound produces the following equation 

Uts — M ts * Cits — Mts Jy-R- Ps,ts 

Many of the IRT test customers are used to working with airspeed in knots, so the following unit 
conversion is needed 


U ts (knots) = 0.592484 • U ts (ft/sec) 


The test section Reynolds number is calculated beginning with determining the air density 


P/5 


^SJS_ 

R • R.s,ts 


The test section air viscosity, (i ts in units of slugs/(ft-sec), is determined from the following equation 


P/5 


Jl.5 

2.270 — 10 -8 

T s ,t 5+198.6 


These parameters are used in the Reynolds number, (10 6 /ft), Re,„ equation seen below 


NAS A/CR— 20 1 0-2 1 6733 


81 



In this equation, the characteristic length, /, is equal to 1 ft. The test section Reynolds number can also be 
calculated based on test article dimension (i.e., chord length) 

Re c 


Table C.l summarizes the inputs, calibrated test section conditions and calculated test section 
conditions as in this report. The program name column lists the variable names for each parameter as they 
appear in the computer subroutine IRTAT. At present, the parameter name and program names do not 
exactly match because the parameter naming nomenclature has been updated to be consistent with naming 
recommendations, detailed in Reference 2. The program names will eventually be updated to be 
consistent with the naming recommendations. Also, note that P air was added as an input to IRTAT as part 
of revision 6. 


TABLE C.l.— SUMMARY OF THE PARAMETERS USED TO DETERMINE 
THE TEST SECTION AERO-THERMAL CONDITIONS IN THE IRT 


Parameter 

Program 

Name 

Description 

Units 

Source 

Average Bellmouth and D-Comer Conditions | 

P T,bm 

POBM 

Bellmouth total pressure 

psia 

Input 

^Pbm 

DPBM 

Bellmouth delta-pressure 

psia 

Input 

Ps,bm 

PSBM 

Bellmouth static pressure 

psia 

Input 

Pl,davg 

TODC 

D-Comer total temperature 

°C 

Input 

p 

1 ratjjm 

PRATBM 

Ratio of bellmouth static to total pressure 

none 

Calculated 

Mbm 

MBM 

Bellmouth Mach number 

None 

Calculated 

P . 

- 1 air 

SBAIRP 

Spraybar air pressure 

psig 

Input 

| Test Section Conditions — Calibrated Parameters | 

Psjs 

PSTS 

Test section static pressure 

psia 

Calibration 

M ls 

MTS 

Test section Mach number 

none 

Calibration 

Tr,ts 

TOTSC 

TOTSF 

TOTS 

Total temperature 

°C 

°F 

R 

Calibration 

(Units) 

(Units) 

Test Section Conditions — Calculated Parameters 1 

Prat 

PRAT 

Ratio of static to total pressure 

none 

Calculated 

Ptjs 

POTS 

Test section total pressure 

psia 

Calculated 

Ts,ts 

TSTS 

Static temperature 

°R 

Calculated 

U ts 

VTS 

VTSK 

Airspeed 

ft/sec 

knots 

Calculated 

(units) 

P ts 

RHOTS 

Air density 

slugs/ft 3 

Calculated 


MUTS 

Air viscosity 

slugs/(ft sec) 

Calculated 


REFT 

Reynolds number 

10 6 /ft 

Calculated 

Re„ 

RECX 

Reynolds number based on model length 

io 5 ! 

Calculated 

y 

GAMMA 

Ratio of specific heats =1.4 

none 

constant 

R 

R 

Gas constant =1716 

lbf-ft/(slug°R) 

Constant 

c x 

CX 

Model chord or other critical dimension 

feet 

PARAM File 
input 
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Appendix D. — 2005 Interim Aero-Thermal Calibration of the IRT 

As mentioned in the Introduction section, there are three basic levels of calibration used in the IRT, 
full, interim and check calibrations. Full and interim calibrations occur in pairs and are designed to 
provide a level of redundancy to ensure consistent operation of the facility. The full and interim 
calibrations occur either following a facility modification or every 5 years if there are no facility 
modifications. The information in this appendix will provide the final test documentation of the 2005 
interim calibration of the IRT including the test matrix, test conditions, data and results. 

An interim calibration uses the same hardware as a full calibration, namely the 9-ft survey rake and 
the 2-D RTD array. For the 2005 interim calibration, testing using the 9-ft survey rake was completed in 
January (January 18, 2005), however, the total temperature calibration measurements using the RTD array 
was not conducted and completed until April (April 1 and April 4, 2005). The time between the 9-ft rake 
and RTD array testing was used for icing cloud calibration and diagnostics. 

The test matrices for the 2005 interim calibration are given in Table D.l and Table D.2. These test 
matrices were designed to: 

• Complete data collection to check the total and static pressure calibration relationships and the 
centerline flow quality using the 9-ft survey rake at the test section centerline at the turntable 
station. 

• Complete data collection to check the total temperature calibration relationship and flow field 
surveys using a two-dimensional array of RTDs. The 2005 test matrix differs from the 2004 full 
calibration by the addition of a spraybar air pressure setting of P air = 40 psig. The addition was 
made to begin a sensitivity study of the spraybar air pressure on the aero-thermal characteristics 
of the test section calibration and flow quality. 

In addition to the fundamental aero-thermal calibration goals, the 2005 interim calibration was used to 
conduct temperature stability studies between the D-comer and test section. The primary concern was the 
overall stability and settling time required to achieve uniform total temperature distribution in D-Corner 
and the test section. The 2-D RTD array was used to monitor the test section temperature distribution in 
conjunction with the D-Corner RTDs. The D-Comer RTDs are facility instrumentation mounted on the 
leading edge of the D-Comer turning vanes. The objective of the surveys was to document the 
settling/stability characteristics to the total temperature profiles in D-Corner and the test section and to 
determine guidelines for IRT operation to insure a uniform and stable total temperature in both D-Comer 
and the test section. The Escort system, the facility data acquisition and display system, was used to 
monitor the standard deviation of the D-Corner total temperature, <5tt,dc> and the time variation in the 
D-Comer average temperature, A T DC . Details and results of the temperature stability study are given in 
(Ref. 3). 


TABLE D.L— TEST MATRIX FOR THE IRT AERO-THERMAL CALIBRATION USING THE 9-ft SURVEY RAKE 
(2005 INTERIM CALIBRATION). 


Rake 
height, Y, 
in. 

Total Temperature, 

A , davgi 

°c 

Test Section Airspeed, U ts , knots 

Spraybar air 
pressure, P air , 
psig 

36 

5 

0, 43, 65, 87, 100, 130, 152, 173, 200, 217, 240, 260, 282, 304, 339 

0 

36 

5 

0, 43, 65, 87, 100, 130, 152, 173, 200, 217, 240, 260, 282, 304, 339 

40 

36 

5 

0, 43, 65, 87, 100, 130, 152, 173, 200, 217, 240, 260, 282, 304, 339 

70 


NAS A/CR— 20 1 0-2 1 6733 


83 




TABLE D.2.— ' TEMPERATURE CALIBRATION TEST MATRIX 
USING THE 2-D RTD ARRAY (2005 INTERIM CALIBRATION). 


Test section 
static 

temperature, 

Ts,ts, 

°C 

Test 

section 

airspeed, 

u ts , 

knots 

D-comer total 
temperature, 

r ^T,davg-> 

°c 

Spraybar 

air 

pressure, 
p . 

1 am 

psig 

Test section 
static 

temperature, 

Ts,ts, 

°C 

Test 

section 

airspeed, 

u ts , 

knots 

D-comer total 
temperature, 

^T.davgi 

°c 

Spraybar 

air 

pressure, 
p . 

1 am 

psig 

-30 

43 

-29.8 

0, 40, 70 

-5.6 

43 

-5.3 

0, 40, 70 

-30 

100 

-28.7 

0, 40, 70 

-5.6 

100 

-4.2 

0, 40, 70 

-30 

130 

-27.8 

0, 40, 70 

-5.6 

130 

-3.3 

0, 40, 70 

-30 

200 

-24.5 

0, 40, 70 

-5.6 

200 

-0.3 

0, 40, 70 

-30 

217 

-23.8 

0, 40, 70 

-5.6 

217 

0.7 

0, 40, 70 

-30 

304 

-17.8 

0, 40, 70 

0 

43 

0.2 

0, 40, 70 

-17.8 

43 

-17.6 

0, 40, 70 

0 

100 

1.3 

0, 40, 70 

-17.8 

100 

-16.4 

0, 40, 70 

0 

130 

2.2 

0, 40, 70 

-17.8 

130 

-15.6 

0, 40, 70 

-0.6 

200 

4.7 

0, 40, 70 

-17.8 

200 

-12.5 

0, 40, 70 

-1.7 

217 

4.6 

0, 40, 70 

-17.8 

217 

-11.6 

0, 40, 70 

-7.2 

304 

4.9 

0, 40, 70 

-17.8 

304 

-5.6 

0, 40, 70 

4.4 

43 

4.7 

0, 40, 70 


3.3 

100 

4.7 

0, 40, 70 


D.l Total Temperature Calibration 

The total temperature calibration relationship is one of a set of three calibrations used to fully define 
the aero-thermal operating conditions in the IRT test section (the other two calibration relationships are 
static pressure and Mach number). The total temperature calibration simply relates the test section total 
temperature measured by the 2-D RTD (average of 49 measurements) array to the total temperature 
measured in D comer (average of 24 measurements). Prior to 2005, a single curve-fit was used to 
represent the total temperature calibration across the operating range of the IRT. This single curve-fit did 
not include data collected at low airspeed with the spray bar air turned on due to the amount of scatter in 
these data relative to the data collected at higher airspeeds. 

During the 2005 interim calibration, an analysis was done to determine if a family of calibration 
curves was needed for the total temperature to account for airspeed and spraybar air pressure effects. The 
final conclusion was that there is a break-point between 50 and 100 knots such that above 100 knots, 
spraybar air pressure no longer has a significant effect on the total temperature calibration. The result was 
to create a set of four calibration curves, one for airspeeds greater than 50 knots and three for airspeed less 
than 50 knots based on the spraybar air pressure setting (0, 40 or 70 psig). The resulting total temperature 
calibration curves are shown in Figure D.l. Table D.3 provides the details for the calculation of the 
calibrated test section total temperature from the measured D-corner total temperature. In order to provide 
some degree of vernier on the spraybar air pressure at the low speed conditions, the slope and intercept 
from each of the low speed calibration curve-fits were plotted against the spraybar air pressure so that 
calibration curve coefficients could be estimated for other spraybar air settings (Figure D.2). The resulting 
curve-fit equations for estimating C T t, o and C T t, i for U ts < 50 knots are listed below. 

For all P a i r settings not specifically listed in Table D.3: 

Ctt,o = 0.06550710 + 0.01720280 • P air + 0.00004323 • P air 2 


For 0 < P air < 40 psig: 


C TT \ = 0.98917044-0.00001881 • P air 
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For 40 < P air < 70 psig: 


C TT l = 0.99428850 - 0.00014676 • P air 

These calculated values for C TT $ and C T t,\ are then used in the calibration equation (Table D.3) to 
determine T Tits at the low speed conditions. 

Based on the 2005 interim calibration, the subroutine used to determine the aero-thermal operating 
conditions in the test section (IRTAT) was updated (revision 6, dated April 11, 2005). The changes were 
the inclusion of the spraybar air pressure, P ain as an input to IRTAT (the computer code name used was 
SBAIRP) and to generalize the equation for calibrated test section total temperature 

Tr,ts ~ Or, o + Qr,i * T T,davg 

Where C T t,o and C T ta are determined as defined above and in Table D.3. 


TABLE D.3. —GENERAL EQUATION AND CURVE FIT COEFFICIENTS USED IN THE 
CALCULATION OF THE CALIBRATED TEST SECTION TOTAL TEMPERATURE. 


Tr,ts ~ Ctt, o + Ctt, i * 7 r,davg 

Test section airspeed, 

Spraybar air 

Curve fit intercept 

Curve fit slope term, 

U ts , 

pressure, 

term, 


knots 

p . 

1 am 

Ctt, o? 

°c/ °c 


psig 

°C 


>50 

>0 

-0.19470324 

0.99014859 

<50 

0 

0.06550710 

0.98917044 

<50 

40 

0.82278309 

0.98841811 

<50 

70 

1.48151795 

0.98401532 


D.2 Flow Quality Comparison 

The interim calibration data are also used to check the consistency in the test section flow quality 
from the previous full aero-thermal calibration. For most of the flow quality parameters this becomes a 
comparison of centerline data since the 9-ft survey rake is only mounted on centerline for the interim 
calibration; however, for total temperature the entire test section cross section is compared as the RTD 
array provides complete coverage during both the interim and full calibrations. 

Figure D.3 shows the centerline Mach number comparison data for the 2004 full aero-thermal 
calibration and the 2005 interim aero-thermal calibration. The Mach number data collected from the 9-ft 
survey rake were normalized by the bellmouth Mach number. Data were also included to show the 
spraybar air injection effects. They-axis scale on the Mach number ratio plots was expanded to show very 
minute detail in the flow field. A AM band was added to each plot to show the absolute variation in the 
data. 

In general, there is no significant difference in the centerline Mach number ratio distributions between 
the 2004 and 2005 calibration tests. For most every case there is nearly an exact match core portion of the 
test section between the data in the sets from each year for the Mach number and spraybar air setting. 
There are two points of slight difference between the 2004 and 2005 data. First is near the inside test 
section wall where the 2004 data is consistently lower than the 2005 data. However, the difference is only 
0.001 to 0.003 in Mach number. The other area is in the right-center of the test section (60 to 90 in. from 
the inside wall) where the 2004 data is slightly higher than the 2005 data, but again by only 0.001 to 
0.003 Mach number. Overall, the 2005 data shows a flatter Mach number profile for all test conditions. 
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Total temperature survey data measured using the 2-D RTD Array is shown in Figure D.4. In general 
the trends in the total temperature data are fairly consistent between the 2004 and 2005 tests. There is 
some improvement seen in the temperature distributions in terms of reduced spread in the data measured 
during the 2005 test, particularly at lower airspeed settings. Since there were not changes to the tunnel, 
this improvement is attributed primarily to allowing more time for the temperature conditions to stabilize. 
Figure D.5 compares the standard deviation in the total temperature measurements over the survey plane. 
This data also indicates little difference between the 2004 and 2005 data sets, except again for the low 
speed settings. 

Figure D.6 shows the resultant flow angle vectors collected using the 9-ft survey rake at 36-in. above 
the test section floor (centerline). This data provides a comparison between the 2004 and 2005 test results 
over the airspeed range of the IRT and at spraybar air pressure settings of 0 an 70 psig; in addition, P air of 
40 psig data collected in 2005 is also included. 

A qualitative inspection of the flow angularity data indicates no significant differences in test section 
flow quality between 2004 and 2005, although the flow angle vector magnitude does appear larger in the 
2005 data. In general, both the magnitude and direction of the flow angularity vectors is consistent 
between the 2004 and 2005 data for like test conditions, particularly for airspeed settings of 100 knots and 
higher. There are some noticeable differences in the flow angle vectors at the lower airspeed settings 
(below 100 knots), but at these low airspeeds, the flow angle probes are not as accurate. 

To provide a more qualitative comparison, average values and standard deviations were computed for 
the pitch and yaw flow angle components (Figure D.7 and Figure D.8). The data were reviewed first by 
using the entire data set (1 1 points across the test section) then by looking only at the core data (the seven 
center probes along the survey line). The following points were drawn from the comparison of the 
average pitch and yaw flow angle data (Figure D.7): 

• The average yaw angle data is repeatable across the IRT airspeed range for both the spraybar air 
off and 70 psig pressure settings. Even the large yaw angles recorded at the low airspeeds with 
the spraybar air pressure at 70 psig show a close match between the 2004 and 2005 tests. 

• The core data indicates that the average yaw angle is about -0.02° for airspeed settings of 
100 knots and higher for both the 2004 and 2005 data. 

• There is a mismatch in the average pitch flow angle magnitude between the 2004 and 205 tests 
seen in both the full-span and core data sets. The core data from 2004 indicate that the average 
pitch flow angle is about 0.1° over the airspeed range; however the 2005 data shows that the 
average pitch angle increases form about 0.1° at 100 knots to about 0.5° at 300 knots. The cause 
for this change in average pitch angle was not determined. 

• Spraybar air effects do not have a significant effect on the average flow angle values (pitch or 
yaw angles) for airspeeds of 130 knots and higher. 

Evaluation of the standard deviation in the flow angle data indicates that there were no significant 
changes in the amount of variation in the pitch and yaw flow angle distribution between 2004 and 2005. 
The pitch and yaw standard deviations from 2004 and 2005 are very consistent over the airspeed range of 
the IRT (Figure D.8). These data also indicate for most the cases; the variation in flow angle is much 
lower in the core of the test section. This is particularly true for the yaw angle data with the spraybar air 
pressure set at 70 psig and the pitch angle data with the spraybar air off. The yaw angle deviation with 
spraybar air off is essentially the same for the full-span and core data. 

Comparison of turbulence measurements made in 2004 and 2005 are shown in Figure D.9. Similar 
turbulence levels and trends over the airspeed range are noted between the two data sets. The only 
discrepancy is at the lowest airspeed setting with spraybar air injection, where the 2005 measurements 
indicated a much higher turbulence level. 
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D.3 Calibration Curve Comparison 

Figure D.10 shows the aero-thermal calibration curves with flow field data from the centerline of the 
test section from both the 2004 and 2005 tests. The calibration curves were developed based on the 2004 
data. The flow field data is the average of the center seven probes on the 9-ft survey rake with the rake 
positioned at the test section centerline (vertical height of 36 in. above the test section floor). There is a 
very favorable comparison for all three calibration curves (Mach number, static pressure and total 
temperature). Figure D. 10(c) also provides a comparison to the total temperature calibration curves based 
on the 2004 and 2005 calibrations. The 2004 calibration curve was based on entire data set (covered the 
entire airspeed range), while the 2005 calibration uses a family of curves to model the total temperature 
and take into account changes in the model over the airspeed range. For simplicity, only the 2005 total 
temperature calibration curve for airspeed greater than 50 knots is included in the comparison. 


o 
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slope and intercept values for spray bar air pressure settings between 0 and 40 and between 40 and 70 psig. 
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Distance from inside tunnel wall, y , inches Distance from inside tunnel wall, y , inches 

Figure D.3. — Continued, (e) M n0 m= 0.200, U n0 m- 130 knots, (f) M n0 m= 0.234, U n0 m- 150 knots, (g) M n0 m= 0.266, U n0 m- 175 knots, (h) M n0 m= 0.306, 
Unom= 200 knots. 
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Figure D.4. — Continued, (h) U ts = 116.46 mph (101.2 knots), T T ,davg=-28A °C, Pair 1 0 psig. (i) U ts = 118.76 mph (103.2 knots), T T ,davg=-28A °C, P ai r7 0 
psig. (j) l/f S =1 16.46 mph (101.2 knots), 7“r,dai/g=-16.3 °C, Pair 1 0 psig. (k) U ts = 118.30 mph (102.8 knots), 7“r,dai/g=-16.3 °C, Pair-7 0 psig. 
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Distance from inside tunnel wall, y, inches Distance from inside tunnel wall, y, inches 

Figure D.4. — Continued, (x) U ts = 253.11 mph (219.95 knots), T T ,davg=0.67 °C, P air = 0 psig. (y) Uts= 292.07 mph (253.8 knots), 
Tr,davg= 0.8 °C, Pair 2 0 psig. (z) l/£s=291.61 mph (253.4 knots), TT,davg=0.7 °C, Pair-7 0 psig. 
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Figure D.5. — IRT test section total temperature standard deviation data collected using the 2-D RTD Array at the center line 
of the RTD Array. 
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Figure D.5. — Concluded. IRT test section total temperature standard deviation data collected using the 2-D RTD 
Array for the full RTD Array. 
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